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ABSTRACT 


The objective of this thesis is to investigate and 
evaluate the effectiveness of modern eStimation methods 
with different array geometries as they apply to the 
problem of bearing estimation. The algorithms were 
selected from those that apply to multidimensional case, 
including MUSIC, PHD, minimum norm, and Capon’s beam- 
former. These four techniques were chosen based on their 
high resolution capability, and their ability to deal with 
eee ee non-uniform arrays and to estimate both 
azimuth and elevation angle of arrival (AOA). Computer 
Simulations were run for linear arrays, circular arrays, 
and combinations of the two. The test conditions included 
(1) two closely spaced emitters, and (2) various levels of 


additive white GausSlan noise. 
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Ae INTRODUCTION 
A. OVERVIEW 


The angle of arrival (AOA) is a valuable parameter to 
be used in deinterleaving signals because a target can not 
change its position rapidly. Even an airbone radar can not 
Significantly change its position in the few milliseconds of 
the waveform pulse repetition interval (PRI) time. As a 
result, the AOA measured by an intercept receiver on the 
radar is relatively stable value. It requires that a large 
number of antenna elements and receivers must be matched, 
elther in amplitude or in phase in order to measure the AOA. 

Direction finding (DF) is the area of electronic- 
warfare support whose objective is to find the AOA [Ref. 1]. 
DF systems pasSively examine the spectrum use by hostile 
emitters and process the signals to obtain enemy AOA 
information. 

The two common approaches to measure AOA are based on 
amplitude and phase comparison [Ref. ile. The phase 
comparison system usually can generate AOA resolution of +1 
degree, which satisfies the modern EW requirement. The 
number of antennas that can be used in an AOA system is very 
low in comparison with electronically scanning antennas in a 
radar system. In an airbone system, the maximum number of 
elements might be ten. In a ship-based system, the number 


might be larger. 


For the phase comparison method at least two separate 
antenna elements with a predetermined space between them are 
required. Since the distances from the emitter to the two 
elements are not the same (except in the broadside case), 
the incident wave arrives at the two elements) after 
traveling uneven paths, and thus it arrives with a different 
phase. The phase delay 1S proportional to the antenna 
spacing, wavelength and angle of incidence. Since the first 
two factors are constant, the direction information can be 
obtained. 

DF systems often receive signals from several sources 
in the same frequency range at the same time and from many 
different angles of arrival. To resolve multicomponent 
wavefields, modern techniques of spectral estimation are 
used. Most recent digital AOA studies have concentrated on 
high-resolution approaches, such as the MUSIC, ESPRIT and 
minimum norm methods ([Refs. 2-4]. In this thesis the 
performance of several modern estimation methods is 
evaluated for different array geometries. This research 
deals with both the one- and two-dimensional problems: 
azimuth angle of arrival for linear arrays, and azimuth and 
elevation angles for a planar array. In addition, an 
underlying assumption is that the receiver is far enough 
from the transmitter and the antenna interactions with the 
platform are controlled so that the arriving wavefront is 


essentially planar. 


The main function of a direction finder is to determine 
the direction of arrival of an incident electromagnetic wave 
relative to the coordinates of the direction finding site 
[Ref. 1]. Figure 1 shows a representative direction finding 
spatial coordinate system with the DF antenna located at the 
origin. The angle of arrival is specified by the azimuth 
angle measured from the x-axis, and the elevation angle 


measured from the x-y plane. Under ideal circumstances the 


electromagnetic field incident on the DF antenna exhibits a 
locally plane wave structure with linear polarization. The 


direction of propagation is indicated by Poynting's 


— 


VeEClor 


P= Ret E\xH, } eu) 


where * denotes complex conjugation, E, is the incident 


electric field intensity, H, is the incident magnetic field 
intensity, 
kxE, 


H, = ' (2) 
Z 





0 


a 


where k is unit vector in the direction of propagation and 


Z 


5 2S JWMpedance Of free space. 


In practice, the incident electromagnetic field is 
usually nonplanar with phase-front distortion created by 


multipath and Scattering, and depolarization produced’ “by 


nonuniform ionospheric propagation effects if such a path 


exists. A full-capacity generic DF should measure angle of 





Figure 1. Direction-finding spatial coordinate system with 
the array located at the origin 


arrival ial three-dimensional space. However, many 
operational situations require only measurement of the 
direction component in the azimuth plane. 

A plane wave, incident at an angle other than the 
normal to the baseline between the two elements, arrives at 
one element earlier than the other. Figure 2 illustrates the 


basic phase delay concept. An incident plane wave arrives at 
an angle @ at antenna 1 inducing a voltage which can be 
expressed aig complex notation as Viexp(jos) . Ritter 


propagrating the distance dsin0, the incident plane wave 


induces a voltage Il’,exp(jor-—@) in antenna 2, where the phase 
delay given by 

®O = (2dn/2)sinO . (2) 
Therefore the bearing angle is encoded as a function of 
phase delay. For the phase difference technique, the bearing 


angle is computed by using Eq. (3), where phase difference 


ra plane 
wave 


Ep front 


- antenna 2 antenna | 
\/ Wy, 


Figure 2. Phase-difference DF Parameters 


is measured and the baseline distance d and wavelength 4 
are known. 

In the case of linear or planar arrays, a phase 
difference can be measured or computed for each combination 
of elements in the array. Then an array correlation matrix 
can be defined and used in one of several spectral 


estimation methods to estimate the AOA. 


5 


B.. THESIS OUTLINE 


Chapter II summarizes several common antenna array 
geometries for one-dimensional DF, including uniform and 
non-uniform linear arrays. The GombinatrionVor cirenlar and 
linear arrays for 2-D DF is also investigated. Chapter III 
covers the .mathematical development of the important 


estimation methods, including multiple signal classification 


(MUS EG) = the estimation of signal parameter via rotational 
invariance technique (ESPRIT), minimum norm, conventional 
and Capon's beamforming, and Pisarenko harmonic 


decomposition (PHD). Also, DF for multiple emitter AOAs is 
introduced by using the array covariance matrix. Chapter IV 
describes the development of a computer code to simulate 
direction finding antenna performance for signals with white 
noise present. Matlab is used to implement the estimation 
methods and perform the array matrix calculations. A 
graphical user interface (GUI) was developed to allow a 
variety of parameters to be changed with a minimum of effort 
on the part of the user. Chapter V shows the simulation 
results, and validates the code by comparing the data with 
previously published results. Finally, Chapter VI presents a 


summary, conclusion and recommendations for future work. 


Ot OS ANTENNA ARRAYS 


Antennas can be viewed as spatial filters; they enhance 
Signals in desired directions while Simultaneously 
suppressing signals in other directions. In a phased array, 
the selectivity 1s achieved by phase shifting and then 
Superimposing the outputs of all elements [Ref. 1]. For a 
passive array, the improvement in the signal-to-noise ratio 
(SNR) 1S potentially equal to the number of elements. In an 
array of identical elements, there are two related factors 
that must be considered. The first is the geometrical 
arrangement of the array (linear, circular, etc.). Second is 
the relative displacement between the elements, once the 
element arrangement is set. 


A. UNIFORM LINEAR ARRAY 


In most EW applications, it is desirable to use as few 
elements as possible. For example, DF arrays on a ship are 
linear arrays placed in the horizontal plane to measure the 
azimuth angle. One linear array. can theoretically cover up 
to 180 degrees of azimuth angle, although it is often 
limited to 120 degrees to avoid operating in the end-fire 
mode [Ref. 1]. I£ the elevation angle is also of interest, a 
linear array in the vertical direction has to be added. As 
far as the AOA measurement iS concerned, the two linear 
arrays are usually treated separately. 


Figure 3 shows a linear array with N elements indexed 


Erom @=U0,1, 2:5, N=1 etong the Xidirection, Mee ies aeoune 


that there is only one plane wave of frequency f/f and that a 


wavefront arrives as shown. If the input signal is a 


Sinusoidal wave, then the qth element will have the output 


2(4,U) = Acos[2nf(t-T, J] 
(4) 


where A is the amplitude of the incident signal, and es 
the time delay at the qth element with respect to the 


reference antenna element (q=0). Since the distance between 
the antenna elements is d, the delay time can be written as 


—qdsin® 
— 


ey) 
C 


where 98 is the incident angle of the input signal, and is 


the speed of light. 





Figure 3. A linear array with a plane wave incident 


The Minus Sign Vin Eo.%5) iso-due to vwthe “fact that. the 
equiphase plane arrives at the gth element before it 
beaches The reterence element, cubvstituting Eq. 915) “inte Eo, 
(4) gives 


qd sin® 
8(q,) = Acos[2fn(t + ————)] . (6) 
C 


~ 


Let kA be a unit vector pointing in the direction of the 


wave propagation and x a unit vector along the direction of 


array axis. Eq. (4) can be written as 


2(q,t) = Acos(2ftn - 2nfgk ae (on) 
G 


~ . TC : 
where k-x= cos(> + 9) =-—sin@ . 


If there are Msignals, the output from the gqth 
antenna element 1S a superposition of all the received 


Signals 


M 2nf.gdk,, x 
ago Ale onpe ae (8) 
Cc 


m=l 


A 


where f, and k, are the frequency and direction of the mth 


Signal. Often exponential form is used to express this 


result as 


2nf_qdk., -% 
ee (9) 


A 
2(4,t) = y A, exp(j(2f ¢t 


m= 


In order to obtain digital data, the output is sampled at 


discrete times and therefore the time ¢ will be replaced by 


integers n=0,/,...,b-1, where b is the number of time 


sample points. If we assume the frequencies of all signals 


2 t 
are the same, the factor e?™™ 


can be suppressed, and a 
Standard phasor quantity results. 

For an AOA measurement, the elements cannot be spaced 
too far apart or an ambiguity will result [Ref. 1]. The 
shortest distance between two antennas must be less than 
half of the wavelength of the highest frequency in order to 
fulfill the Nyquist theorem [Ref. 1]. Therefore, one must 


sample the incident wavefront twice per cycle. 


5B. NON-UNIFORM LINEAR ARRAY (MINIMUM-REDUNDANCY) 


By making use of a theorem by Caratheodory [Ref. 5] it 


is shown that, ‘for a given number of elements WN, there 
exits a distribution of element positions which results in 
superior spatial spectrum estimators than the uniform linear 
array. In addition to less hardware, certain nonuniform 
arrays are capable of large dynamic range (spectral peak to 
background level), lower sidelobes, and the relatively small 
estimate bias values. A special type of nonuniform array is 
the so-called minimum redundancy array (MRA), in which 


integer spacings of the base spacing occur only once. A MRA 


is capable of resolving MWN-1/)//2 multiple (simultaneous) 


signals, whereas a conventional array can only resolve N-/. 
Some minimum-redundancy array configurations are shown in 


Table 1. A dot in column 3 represents an element and the 
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integers the number of base spacings between them. (The base 
spacing is typically A/2.) 

A disadvantage of minimum-redundancy array is that its 
resolution is not easily increased except by adding more 
antennas and rearranging the array to the optimum 


configuration for the new number of antennas. If there are 


M signals, the output from the gth antenna element is 


M Mf nd ky * X 
2(4,)) = x Am &XPG(20f,,¢ - ————_)) (10) 


where d, 2S the Cistance between first antenna ance he wou 


antenna. 


Table 1. Some minimum redundancy array configurations. 





©: CIRCULAR ARRAY 


A planar array can be used to measure both the azimuth 
and elevation angles. These types of arrays occupy more 
space than the linear array but they might be necessary for 
some special aircraft and shipboard application. Planar 
arrays include circular and rectangular (of which crossed 


linear arrays can be considered a subset). A circular 


deal 


antenna array is shown in Figure 4. The radius of the circle 
is Rk, and there are N elements from qg=0 to N-I/. The 
array 1s in the x-y plane and the first antenna element is 
the reference element (g=0) and it is on the x-axis. The 
location of the elements are at integer multiples of the 
angle 6,=2n/N. 

Assume that a plane wave with amplitude 4, and 
frequency f/f, is incident on this array with azimuth angle 9, 
and 9, (the angle between the incident ray and the normal of 
the array). The output induced at qth element can be 


written as 


9(q,t) = A, exp(i(2nf,(t- yy (12) 





Figure 4. Circular antenna array 


where the space delay T referenced to the center of the 


ql 
array is 


eZ 


t,, =—K, sin9, cosg,, . em) 
In this case, because the antenna locations are equally 


spaced in azimuth, the angle 0,, can be written as 


ql 


941 = 99, -9; . ie) 


Therefore 


: 21q 
K, sind, Co —,) 


8(q,0) = A, exp(J2mf (t+ )) 


, 21g 
Rk,sin8, ar ~,) 


= A,exp(j2n(f,t + 7 Ds (14) 
/ 


The incident direction of the input signal can be 


A 


represented by a unit vector k, in Cartesian coordinates as 
k, =—(sin@, coso x +sin@, sing, y+cos8,2 (15) 
where xX, py and = are unit vectors in the Cartesian system. 


The minus sign in front means the vector is pointing to the 


©rigin. The postion vector GCE the Quill antenna yelement. 16 


r. where 


r, =cos(qd, )x +sin(qd, )y (16) 


so that 
k,-F, =—sin®, cos(qd, —)) . (ys) 


Extending this to cover multiple signals (m=I,2..,.M) gives 


abe: 


aA 
am 


M . eee 
2(g,t) = > A, OU army é (18) 


m=] ae 
The linear array can only provide azimuth angle 
information because the corresponding propagation vector is 
two-dimensional. The circular array can provide both azimuth 
and elevation information since the GOnneSspoOnaing 


propagation vector is three-dimensional. 


D. SPIRAL ARRAY 


We examine a variation of the circular array, Eq. (18), 


where the radial distances of the elements increase in steps 
Ol ere) ees eae Bense e yg ) eeee eee Clea = EO omeem 
qth antenna element and the center point. The element 
locations lie on an Archemedian spiral. In a circular array, 
aii anbenna (elements fave Ene (same =i - The “outpur "or aan 


antenna element is as follows: 


Rk.-? 


M 
2(g,t)= XA, expG2n(f,,¢--+* + )) . (19) 
m=] dn 
E. COMBINATION OF LINEAR AND CIRCULAR ARRAYS 


A three dimensional array can be constructed by 
combining planar and linear arrays. Figure 5 shows a 


circular array in the x-y plane and a linear array along 


the =z axis. We examine the special case where the elements 
of the circular array have 6,= 90°, which results in two 
elements on the x and y axes equally spaced d from the 
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center. Let there be 3 the elements along the z axis, also 


spaced d. Therefore, the array has seven elements, all 


equally spaced. 





x 


Figure 5. There dimensional seven element array 


The positions of the 3 elements of linear array are 


given by the position vector 
r, =(h-l)dz 
where h=1,2,5, so that 
kT, =—c0s(0, (h—-L)d) . (20) 


If the array receives M signals, the output of Ath 


vertical antenna element of the linear array is 


Binet ) = _ x A, SEXPU(2 1, Ls 2a th a: ) CZ ay) 


“fi 


Ra) 


Another four antenna elements comprise the planar array 


which has an output of the form 


A 


M . hei 
§ circular (q, t) — _ ae exp 2n(f,,f ~ y )) (22 ) 
where qge=/,...,4, M is number of signals, and k,-F is the 
same as Eq. (18). This array combines linear and circular 
arrays to provide 2-D direction finding. From Eq. (21), it 


is evident that the vertical array only provides elevation 
angle information. 

We have described the various array geometries. The 
elements comprise the data acquisition system which samples 
the incident signal wavefronts. The next step is to carry 
out the estimation process. In the next chapter, we 
introduce spectral-based algorithmic solutions to estimate 


the AOA. 
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Ei. SPATIAL SPECTRUM ESTIMATION 
Many different high-resolution approaches can be used 
to estimate AOAs from digitized input data. In this chapter, 
five high resolution methods will be discussed: 
(1) Multiple Signal Classification (MUSIC) [Ref. 2]. 
(2) Estimation of Signal Parameter via Rotational 
invariance (ESPRIT) [Ref. 3]. 


(3) Capon's beamformer (Maximum Likelihood Method) 


[Ref. 6]. 
(4) Pisarenko Harmonic Decomposition (PHD) [Ref. 7]. 
Co Minimum norm [Ref.1]. 
All of these are cabable of (1) high resolution, 2) 
estimating both azimuth and elevation AOA, and (3) handling 


three-dimensional non-uniformly spaced array signals. 


A. MUSIC 


MUSIC is a technique used to determine the parameters 
of multiple wavefronts arriving at an antenna array from 
measurements made on the signals received at the individual 
elements [Ref. 2]. 

7 The waveforms received at the N array elements are 
linear combinations of the incident wavefronts from M 
narrowband signals and noise. Thus, the multiple signal 


classification approach begins with the following model for 


characterizing the receive vector G as 


dy 


GC=aAray (23, 


where V ais noise vector. Expanding Eq. (23) 
§] F, V, 
E2 aj(9;) a)(9z) ... @)(Ox,) fy V; 


| 42(9;) 4,(87) (yy) 


fs Need Ue 
EN _ Fuel |Yn 


The size of each matrix is as Follows: 





Ge Veceor jon 2 ememe -ouroutSs . wy oo) | 

A: matrix of mode vectors, N by M 

i 2 TeeChOr O8 St Gna lsat aon = | 

VY: noise vector, N by 1 

The incident signals are AG Oe need in amplitude and 

phase at some arbitrary reference point, for instance, the 
Origin of the coordinate system, by the complex VeECEO Ls 
The elements of G and A are also complex in general. The 


a,{8,) are functions of the signal arrival angles and the 


array element locations. That is, a,(®@,) depends on the ath 


array element, its position relative to the origin of the 
coordinate system, and its response to a signal incident 
Prom the direction of the mth source. The mth «column “of 4 


can be considered as a ‘mode’ vector of responses to the 


Gibeection of arrival 0, of the mth-signal. This WV by 1 


mode vector will be denoted by a(@). 
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The M by M covariance matrix of the G vector is 


4 Seri eer a * St aia 
S=E[GG |=GG =AFF A +VV (25) 
where E is expectation operator and the asterisk denotes 


complex conjugation. ce (29) can be rewritten as 


— 
S=APA +AS,, where P=FF is a diagonal M by M matrix and 
S, is the noise covariance matrix. 


When the number of incident wavefronts M is less than 


the number of array elements N, then APA is singular. It 
has a rank less than N; therefore |APA’| =|S- aS,| = 0 ; Tas 


equation is only satisfied with 2 (not to be confused with 
wavelength) equal to one of the eigenvalues of S in the 


matrix of S&. Therefore A can only be the minimum 


o 
eigenvalue 4,... Observe that any vector orthogonal to A is 
an eigenvector of S with the eigenvalue ao, the noise 
variance. Hence, we can write 

SeURG =U AeG 0 oh cre (26) 
where A=diag{2,,A,,...,44} 1s a diagonal matrix of real 
eigenvalues. The superscript H represents the Hermitian, 
U, is the signal eigenvectors and U, is the noise 


eigenvectors. The noise eigenvectors can be used to form an 


estimator for the spatial spectrum 


nl, 


l 


= (29) 
a(0)" UUs a(8) 


Puy (0) a 


Although MUSIC is highly robust 1t requires echaracreri zation 
of the array response and a Spatial search through all 
possible angles of arrival,@. 

Be ESPRIT 

ESPRIT [Ref. 3] eliminates the undesirable features of 

MUSIC, but only under certain constraints. Assume that there 
are N points of data (N antennas), which can be divided 
into two groups, each consisting of N-/ points. The first 
group has an output signal x() =As(t)+n,7) and the second 
group y(tH=ABs(t)+n,(t), where A is a steering vector (not 
pee EO <A ans Boe (24)), s(t) denotes the baseband 


waveforms, and B is a diagonal matrix 


v, O 0 
0 v, 0 0 
B= 0 (28) 


10 000 wy, | 


my ' 
where v, = exp(-j—dsin8,) and d is antenna spacing. These data 


will be used to make two matrices by the covariance 


approach. 
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Let Rx, be the auto-covariance matrix of x(t and Ryy 


the cross-covariance matrix of x(J) and yd). If there are M 


Signals then 


Rux= E[x(t)x ()]= ARA +07] (29) 


Ryy= Elx(n)y' (Ol= ARB A od, (30) 


* 
where R=(s(t)s (t)] is the signal covariance matrix under the 
assumption of uncorrelated noise, / is the identity matrix 


and J, is an N by N matrix with ones along the first lower 


diagonal off the major diagonal and zeros elsewhere. For a 
pairwise matched array with doublets the cross covariance 


matrix becomes 
Re Ake A (ay 


In the absence of coherent signals, eliminating ore we 


obtain 

Cyx=Rex- ol = ARA 

Cxy=Ruy- 07d, = ARB A’ 
which gives 

Cyx- ¥Cxy= AR(I-yB)A. (22) 
The singular values of Hg@.({32) are given by Ene rooms of 


I —yB"|=0. The desired singular values are 


Ze 


y=VM, , Kel,2,.-,M. 


To find the angles we note that for a uniform linear 
array, the angles of the eigenvalues are equal to 2drsin@ so 


that 


- g 
6, =sin iC) (33) 


— 


where 9 is obtained from y, =exp(j3,). Thus, the direction | 


GF Jarrival can be: Obtained without ~a. -searen technigue. in 


this respect computation and storage costs are reduced 


Cone rderalin. 
ce BEAMFORMING 
ae Conventional Beamformer 
The first attempt EO- automatically @lCcalize So rem 


sources uSing antenna arrays was through beamforming 
technigues. The idea is to steer the array beam in one 
Gdisection, at a timesgeanc. measure mine Woueru ey cower ay! ic 
steering locations which result in maximum power yield the 
DOA estimates. The array response is steered by forming a 


linear combination of the antenna outputs 
Noe ih 
ee Ce) (34) 
i=] 


where gf) is the output of the ith element, and w,(f) the 


corresponding complex weight factor. The output power for }b 


time samples is 


Ze 


ae 2 ] bel 
P(W) = YIyo| = YW9GaG" (NW =Ww "RW (35) 
f=0 f= 


where K is the signal covariance matrix defined previously. 

The conventional beamformer is a natural extension of 
classical Fourier-based spectral analysis to antenna array 
processing. Suppose we wish to maximize the output power 
from a certain direction 90, given a signal s(t) emanating 
from direction @. A measurement at the element output is 
corrupted by additive noise nt) and written as 
g(q=ales(jtni). The problem of Maximezing = Ene ourouL 
power is formulated as, 


max[ EW" G(yG" (NW } J= max[W" EYG(NG" (OWT 7 


= max} E)s(o) W"a(®)| 4 ot) | (36) 


where the assumption of spatially white noise is used. The 


resulting solution is 


a(Q@) . (37) 


Ja" @)a(6) 


The above weight vector can be interpreted as a spatial 


W, 


filter, which has been matched to the impinging signal. 
Using the weighting vector, the classical spatial 


spectrum estimator 1s 
P,- (8) = a” (8)Ra(6) (38) 


where R is the signal covariance matrix. 
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For a uniform linear array of isotropic sensors, the 


steering vector 1s 


a(O0)=|/ Be a » calle ost -10]" (39) 


where @®M=-kdsin0, N is number of elements and T denotes 


transpose. The standard beamwidth [Ref. 6] for a uniform 
. 2 
linear array 1s We! SO sources whose angles are closer than 


this beamwidth will not be resolved by the conventional 
beamformer, regardless of the available data quality. 

Ze Capon’s Beamformer 

In an attempt to alleviate the limitations of the 
conventional beamformer, such as its resolving power of two 
sources spaced closer than a beamwidth, an improved method 


was proposed by Capon [Ref. 7]. The optimization problem was 


posed as min{P(W)}, subject to W"a(8)=1. Hence, Capon’s 
beamformer (also known as the minimum variance 
distorsionless response filter) attempts to minimize the 
power contributed by noise and any signals coming from 
directions other than 90, while maintaining a fixed gain in 
the look direction G77 Ihe opeimal j7 “can be found using 
methods such as the technique of Lagrange multipliers [Ref. 


6], resulting in 


R7'a(6) 


a 40 
a” (8)R~‘a(8) ate 


CAP 
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Using the above weight vector, we can find the following 


Spatial spectrum 


| 


le 0) = —________ 
cap (8) a” (8)R-"a(6) 


ae 


The power minimization can also be interpreted as 
Sacrificing some noise suppression capability for more 
focused rejection in the directions where there are other 
sources present. The spectral leakage from closely spaced 
sources is therefore reduced, though the resolution 
capability of the seen beamformer is still dependent upon 
the array aperture and on the signal-to-noise ratio though 
i 
om PHD 

The PHD is based on the use of a noise subspace eigen- 
vector to estimate the incident angle and assumes that the 
process consists of M complex signals in additive complex 
white noise. It derives the angle of the signals, their 
powers) and the white noise variance from the known 
autocorrelation sequence. This algorithm accounts for the 
conjugate symmetry of the eigenvector. This method is based 
on an estimation of angles by using the orthogonality 
property between the signal vectors and vectors in the noise 
subspace. The minimum eigenvector is orthogonal to the 


Signal vectors, thus we have 
a," (®)uy =0, i=1,2,..M. (42) 
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Onee the Minimum eigenvector yy 15 found the from (7 ame eee 
(26), the angles of the signals may be determined by 
factoring the eigenfilter polynomial [Ref. 7]. The roots are 


guaranteed to be on the unit circle. The a(9) is a steering 


vector that will be zero whenever 9=98., one of the angles 
of the signals. This leads to the angle estimator function 


la” (8 )uy 





Le MINIMUM NORM METHOD 

This method and the MUSIC method are somewhat alike. 
The basic idea is to find a vector D that is a linear 
combination of eigenvectors in the noise subspace. This 


vector D can be written as 

DS, eye ee ees |e Se Opa = Onde (44) 
where the superscript T signifies the transpose and N the 
number of antenna elements. In this equation 6, is assigned 


to be unity. The square of the norm of Eq. (44) is 
2 pa 
ey = 2 5 
i=0 


where M is the number of sources. The minimum norm method 
minimizes this quantity and hence the name. The procedure is 


as follows: 


Zo 


1. Find the eigenvectors U of the covariance matrix R. 
Large eigenvalues correspond to Signals and small 
eigenvalues correspond to noise. 

2. The vector D can be found from the signal subspace U, 
ii BOs (26) 4 

3. The vector D can also be obtained from the noise 


subspace. The signal subspace U, can be partitioned as 


Uno Yor - Uon-] 
H 
u u ou 
10 a Ne i 
Us = = -| : (45) 
U 
Uyo Up - Uyn-] 


where n” is the first row and U, consists of the remaining 


Fows of U.. Inhese two Matrices canes: woleeewuas 


nN” =([Hop Yop + Mow-r] 
ip pe 

ve (46) 
tUyo Uy »- +» + Uyn-7 | 


The D vector can be obtained as 


] 
be “Us (47) 
(1-7n) 


4. Once the vector D is obtained, the angle estimator 


function can be written as 


] 
Pen 0\ eee (48) 
un (@) a(Q)DD" a(6)" 


Zul 


So far several high resolution technigues have been 
applied -to the problem of resolving the directions of 
arrival of incoming Signals. FKOm a °Sstartstical point oe 
view, the signals arriving at the array can be regarded as 
random, and thus the signals and noise are uncorrelated. In 
the next chapter, we implement these methods based on finite 


observations. 
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IV. IMPLEMENTATION AND COMPUTER CODES 
A. PARAMETRIC DATA MODEL AND FLOWCHART 


The AOA estimation methods described in Chapter III 
were implemented in a computer code. Output data were 
gathered for a variety of arrays and signals, and a 
performance comparison made. This chapter describes the 
structure and operation of the computer simulation, which is 
written in Matlab. The approach is model-based in the sense 
that it relies on certain assumptions made on the 
observation data. Since the natural causes responsible for 
Signals and noise are often unrelated, it 1s customary to 
assume that the signals and noise are uncorrelated with each 
other. Figure 6 shows the flowchart of the computer program, 
with the main steps described in the following sections. 


1. Input Parameters 


For an antenna array with MN elements, the size of 


correlation matrix (correlation between ith and jth 


element) is N by N. In this case, the maximum number of 


sources that can be detected cannot exceed (N-1). Because 
the effectiveness of these spatial spectrum methods are a 
function of signal-to-noise (SNR) ratio, we evaluate the 
performances of each antenna array by changing SNR at the 
input. Also, we must consider the spatial resolution 
capability of each estimation method when applied to each 
array geometry. For example, test signals that have the same 


azimuth or elevation angle are generally harder to resolve 


Zo 





1. Number of elements 4. Angles of arrival 
2. Numbcr of sources 5. Element spacing 
3. Sample points 6. SNR 












1.Signal model: Laplacian 
distribution random processing 
2.Complex white noise model 









Uniform Non- 


Linear Uniform 
Array Linear 


Array output matrix 


J 
Singular Value Decomposition 
U 


Spatial spectrum method 
















eR 


ee 


Plot: Angle of arrival vs. Power 


Figure 6. Computer simulation flowchart 


than those not located in a common principal plane. 


examine the effects of the number sampling points. 
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MUSIC PHD Capon’s MLM ESPRIT 










We also 


The more 


sampling points, the more accurate the AOA estimate but it 


occurs at the expense of increased processing time. 


Generally the antenna spacing must be less than 0.54 in 
order to avoid the ambiguity problem. 


2. Signal and Noise Models 


In practice, the white noise and signals are modeled as 
a independent and identically distributed random processes. 
That means we interpret the received signals as the sample 
of some waveform at certain specified instants of time, and 
is referred to as a realization of a discrete time random 
process. In all cases the signals that arrive can be 
regarded as random, we assume that signals sS(m) are 
Laplacian probability distributed processes. Because all 
arrays are not diverse in polarization (1.e., they have all 
elements identically polarized), we assume that the 
frequency and polarization of all received signals are the 
same. Multi-path and re-radiators in the vicinity of 
receiving site which may have a mutual coherence are not 
considered. 

The signal and white noise generation process are 
described by the following Matlab pseudo-code: 

u(n) = rand(n)—0.5 


s(n) = —(u(n)) * log(1 — 2* | u(n) |) * A 


win) = “Evn) 


v(n) = randn(n) + 7 * randn(7) (49) 


emi 


where n=0,1,...,4095, A is the signal amplitude, ie the 


peak amplitude of white noise and the noise wn) is a 


complex Gaussian random variable with both the real and 


imaginary parts having a variance equal to 1. The Matlab 
functions rand and randn can be used to generate Gaussian 
distributed random values. The voltage ratio re NE 


deberminicc= tie QSNR th de ee Mice 2007 eres 


COnvVentence "we lets A =i Ssorvenat 


N,=10"°". (50) 


3. Antenna Array Model 


Based on the discussion in Chapter II the steering 
matrices for five array geometries can be determined. The 


steering matrix for the uniform linear array takes the form 


] Se ] 
A(o) = (51) 


exp(j2m(N -/)dsing,;) + -: : exp(j2n(N - /)d sind ,,) | 
where N is the number of antenna elements, M is the 


—_ 


number of sources and 6, is the AOA of the Mth source, 


SU yi A 


The steering matrix for a the non-uniform linear array 


1s, 


BZ 


A(o) = : . . (52) 


| exp(j2ndy_,sin,) eee ee exp( j2md y_; Sin 4, ) | 
where d, is the spacing between the reference element and 


the (f+1)th antenna. For a minimum redundancy array, the d, 


ave speciti1ed i Table... 


The steering matrix Lor circular array is 
exp(j2mR, cos(—o,)sin0,;) - - + exp(j2nR, cos(—o,)sin9,, ) 

A(0,0) = (53) 
_exp( J2TR, Cos(dy,)sin®,) - - + exp(j27R, cos(d yy, )sin 9 y, )| 

where 9,,=no,-90,, 9, is the elevation angle of mth source, 

>, is the azimuth angle of mth source and k, is the 

radius of the circle. For the present analysis it is 


assumed that Rk, is equal to 0.54, 0S, <2n, and Ces WF, ; 
The steering matrix for the spiral array is 
exp(j27R, cos(—o,)sin8,) - - + exp(j2nk, cos(—o,)sinG,, ) 


A(9,0) = (54) 


Lexp(j2mRy cos(>y,;)sin®;) - + + exp(j2mRy cos(d yyy) sin Oy, ) | 
where R, is the distance between mth antenna element and 


center point. In our model, the maximum R, is 0.54. 


n 


BP 


The steering matrix for the combination linear and 
circular array can be partitioned into £sub-matrices 


associated with the linear and circular parts of the array 


- - ae | 


exp( /27d sin 0,) exp( j27d'sin8,, ) 
exp( /4md sin0,) exp( /4md sin@ 
A(8,) = | ’ ae m) (eS) 
exp(j2md cos(—,)sin@,) exp(/2md cos(—o,, )sin@,, ) 
| | | 
exp(j2mdcos(oy,)sin®;) + - + + exp(j2mdcos(d yay) SiN By, ) | 
where d is 0.4/4. The first three rows in Eq. (55) are for 
the linear array elements. For these terms there is no 


dependence on arrival azimuth angle. 


4. Correlation Matrix and Spatial Spectrum Methods 


The output of the array can be calculated by using the 
Signal, noise, and antenna models in Eq. (24). First the 
Matlab function is used to find the singular value 


decomposition. The result 1s substituted into the various 


algorithms (Eq. (272 (41), (43), or (48)) with the 
appropriate antenna array steering matrix (Eq. (51), (52), 
(53), (54), or (55)). The outputs for MUSIC, PHD, minimum 


norm and Capon's beamformer are angle of arrival vs. power. 
The power must be computed at every angle. We must steer the 
antenna array in all azimuth directions (0 to 360 degrees) 
and elevation direction (0 to 90 degrees) in order to find 


the AOAs. The output of ESPRIT is the direction of arrival 


34 


which can be obtained without a search. This method has only 
been applied to a uniform linear array. 


Be THE IMPLEMENTATION OF ARRAY MATRIX DECOMPOSITION 


The signal measured at the output of any element in the 
array differs from the signal actually received by an amount 
attributed to noise. Thus, we have an observed signal, which 
for each element of the array, consists of the actual 
received signal plus narrow-band noise. Calculations used 
56=4096 time sample points for each antenna element, so the 
size of received signal matrix is N by 4096, where N is 
the number of elements. The size of signal-in-space matrix 
is M by 4096, where M is the number of sources. The 
observed signal matrix is equal to the received signal 
Matrix times the propagation delay matrix. The observed 
signal matrix is the antenna output matrix G, from which we 
find the covariance matrix of the observed signal matrix 
defined by 

Ree 2G Ge 

The covariance matrix R has two important special 
properties. First, it 1s Hermitian, that is, it is equal to 
the conjugate transpose of itself. Second, it is positive 
semi-definite if element noise is present. Thus, for 
computer coding purposes 


R-= £1G Gl 2G Ga. (56) 


fakes) 


The size of R matrix is N by WN. The Matlab function SVD 


1s used for singular value decomposition. It has the form 


[U. S, U]=SVD(R) 


where R=USU" (Eq. (26)). The eigenvector U is an N by N 
orthogonal signal and noise subspace matrix, and S is an 
N by N diagonal matrix. The singular value decomposition 


igs accurately computed by this function. The columns of U 


are the eigenvectors of RR’, and the entries along the 


diagonal of S are the correspondingly ordered nonnegative 
square roots of the eigenvalues of R'R, which are also the 


eigenvalues of RR’. The numbers in S are the Singular 
values of the matrix kK, which is determined by the number 


of sources. After calculation of U and S, they are 
substituted into Ene Spectral estimaroro wy mec med 
(43), and (48) to determine the AOAs. 

The MUSIC and PHD methods make use of the noise 
subspace matrix; the minimum norm method uses the signal 
subspace matrix. Capon's beamformer does not use the 
Singular value decomposition of the covariance matrix. For 
2-D direction finding, it takes a significant amount of CPU 
time to search space in order to find the AOAs. The default 


search range in azimuth angle is 90 degrees. 
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Cr GRAPHICAL USER INTERFACE 


To compare the relative performance of the various 
arrays and algorithms, the simulation must be run for many 
cases. A graphical user interface (GUI) was developed to 
vary the input parameters, as shown in Figure 7. Also 
included is an on-line help feature. The help button 
provides brief explanations of the input parameters. The GUI 
provides default values to demonstrate the angular 
resolution of MUSIC for a uniform linear array and a minimum 


redundancy array with different signal-to-noise ratios. 
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source number: “48 
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Figure 7. Graphic user interface 









In the next chapter, we use this GUI to present a 
performance comparison of superresolution spectral 


a7 


algorithms for elevating and analyzing their effectiveness 


in several antenna and emitter geometries. 
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Vi. SIMULATION RESULTS AND ANALYSIS 


In this chapter the results of the computer simulation 
are presented. The performances of the various algorithms 
with different array and emitter geometries are examined. 
Computer simulations were run for the conditions of: (1) two 
closely spaced AOAs for testing azimuth and elevation 
angular resolution, (2) various levels of additive white 
noise, (3) multiple signal sources, and (4) changes in the 
antenna element spacing. 

By comparing the test results, we can determine which 
combination of array and algorithm operates best under the 
given set of conditions. The power spectra are plotted as 
surfaces. For all cases, the z-axis represents magnitude of 
the power, the x-axis represents as azimuth angle and the 


y-axis the elevation angle. The majority of the cases use 7 


antenna elements; the exceptions are Figures 35 and 37 
which are for 10 elements. 


A. ANGULAR RESOLUTION OF ONE-DIMENSIONAL DF 


nae el simulations were run to determine the 
comparative performance of algorithms in resolving closely 
spaced AOAs. The signals were assumed to be of equal 
frequency, polarization and mutually incoherent. 

Simulations were run for two signals with angles of 
arrival of 15 and 18 degrees; the corresponding spectral 
estimates are shown in Figure 8. The rows correspond to 
spectra estimates for the PHD, MUSIC, minimum norm and 


Capon's beamformer algorithms, and the columns correspond to 
ee 


the uniform array (left) and non-uniform linear array 
(right). From this figure, we can see that the non-uniform 
array can achieve better angular resolution. All algorithms 
except the Capon's beamformer provide sharp distinguishable 
peaks, but MUSIC has a smoother spectrum (the separation of 
the peaks is not aS pronounced). Figure 9 shows” an 
interference simulation for 6 sources. The non-uniform 
(minimum redundancy) array provides better dynamic range for 


all algorithms. 


1S TWO-DIMENSIONAL DIRECTION FINDING 
1. Azimuth Angular Resolution 


A range of SNRS were simulated. Figures 10, 11, and 12 


show the results for a high SNR (=25 dB), at an elevation 


angle of 15 degrees and AAZ=5 degrees. The MUSIC, PHD and 
min-norm methods have similar results for all three 2-D 
array geometries axeene the Capon's beamformer. However, 
when the SNR is lower (15 a@B), the circular array performs 
better, aS indicated by a comparison of Figures 13, 14, and 
15. Even though two peaks exist, there is no notch between 
them and therefore it may be difficult to correctly resolve 
the two spectral peaks. Figures 16, 17, and 18 show the 


results for SNR=30 daB, and two sources at the same elevation 


angle but separated in azimuth (AAZ=3 degrees). All arrays 
can resolve the two peaks. However, for the lower SNR level 


case, the circular array 1S generally superior. 
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Figure 8. Comparison of uniform (left) and non-uniform (right) 
linear arrays, 2 sources, 7 elements, separation angle 3 degrees 
(15° and 18°). SNR=25 dB 
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Figure 9. Unifrom (left) and non-uniform (right) linear arrays, 
SNR=25 dB, 7 antennas, 6 sources: 15°, 25°,35°, 45°, 55°, and 65° 
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Figure 12. Spiral array, SNR=25 dB, 2 sources at (16°,15°) and 
(21°,15°), 7 elements 
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Figure 13. Circular array, SNR=15 dB, 2 sources at (16°,15°) 
and (21°,15°), 7 elements 
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Figure 14. Lineart+circular array, SNR=15 dB, 2 sources at 
(16°,15°) and (21°,16°), 7 elements 


MUSIC PHD 


power(dB) 
NO 
© 
power(dB) 






a 


Bie ae : 50 
EL(deg) O 0 ~AZ(deg) El(deg) 9 O AZ(deg) 
Mini-norm Capon MLM 





power(dB) 





400 


50° ~“E0 
EL(deg) O 0 AZ(deg) El(deg) OU 0 AZ(deg) 


Figure 15. Spiral array, SNR=15 dB, 2 sources at (16°,15°) 
and (21°,15°), 7 elements 
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Figure 16. Circular array, SNR=30 dB, 2 sources at (16°,15°) and 
(19°,15°), 7 elements 
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Figure 17. Lineart+circular array, SNR=30 dB, 2 sources at 
(16°,15°) and (19°,15°). 7 elements 
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Figure 18. Spiral array, SNR=30, 2 sources at (16°,15°) and 
(19°,15°), 7 elements 
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Figure 19. Circular array, SNR=25, 2 sources at (16°,15°) and 
(16°,18°), 7 elements 
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2. Elevation Angular Resolution 


Simulations were run for the two signal case with 
angles of arrival at 18 and 15 degrees in elevation and 16 
degrees in azimuth. For the SNR equal to 25 dB, the 
corresponding spectral estimates are shown in Figures 19, 
20, and 21. The results for the circular and linear+circular 
arrays are almost the same, but the spiral array cannot 
correctly resolves the peaks. For a 15 dB SNR the data, 
which is shown in Figures 22 and 23, indicate that all 
algorithms, except Capon's beamformer, correctly resolve the 
two spectral peaks with a circular array. Only two 
algorithms can resolve the peaks for the circular array, so 
the elevation angular resolution of the linear+circular 
array appears to be better for this scenario. 


3. Multiple Signal Sources 


Multiple signal simulations are shown in Figures 24, 
25, and 26. There are four signals arriving at azimuth and 
elévationm angles of (10°,20°), (55°, 1547 (30 Alc and (457 ,-25°) 
with equal amplitudes. The circular antenna array and 
linear+circular array can resolve the four peaks but the 
Circular array has more sidelobes. Only the MUSIC algorithm 
correctly resolves the signal sources for the spiral array. 


In this scenario the linear+circular array is best. 
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Figure 20. Linear+circular array, 
(16°,15°) and (16°,18°), 7 elements 
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Figure 21. Spiral array, SNR=25 dB, 2 sources at (16°,15°) 


and (16°,18°), 7 elements 
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Figure 22. Circular array, SNR=15 dB, 2 sources at (16°,15°) 
and (16°,18°), 7 elements 
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Figure 23. Lineart+tcircular array, SNR=15 dB, 2 sources at 
(16°,15°) and (16°,18°), 7 elements 
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Figure 24. Lineart+circular array, SNR=25 dB, 7 elements, 4 
sources at (10°,20°) , (55°,15°) , (30°,60°), and (45°,25°) 
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Figure 25. Circular array, SNR=25 dB, 7 elements, 4 sources at 
(10° 720 0 (55 715 ens 0 260 ) ana (45 425 ) 
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Figure 26. Spiral circular array, SNR=25 dB, 7 elements, 4 
sources at (10°,20°), (55°,15°) , (30°,60°), and (45°,25°) 
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Figure 27. Circular array, SNR=15 dB, 7 elements, 2 sources 
at (16°,15°), and (16°,20°) 
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4. Comparison of Azimuth and Elevation Angular 
Resolution 


The estimation methods and arrays were compared for 


two emitters with a 5 degree separation. Figures 27, 28, and 
29 are for SNR=15 dB and AEL=5 degrees; Figures 13, 14, and 


15 for SNR=25 dB, AAZ=5 degrees. We see that the elevation 
angular resolution is better than the azimuth angular 
resolution since the sidelobes in azimuth are higher than 
the sidelobes ima, elevation ea all three array 
configurations. Even at lower SNRs, all three arrays 
employing the MUSIC method can distinguish the two sources 


in elevation. 


5. Improvement of Angular Resolution of Two-dimensional 
DF 


Next we investigate changing the radius of the circular 
array to 1.5 wavelengths. The results are shown in Figures 
30 and 31. Comparing Figures 13 and 14 with Figures 30 and 
31 (SNR=15 dB, AAZ=5 degrees), the angular resolution is 
better with the larger radius as eee Larger spacing 
may be susceptible to ambiguities. From Figures 32 and 33, 
where the range of azimuth angles is from 0 to 360 degrees, 
it is evident that there is no serious ambiguity problem. 
Thus, we can improve the angular resolution by.changing the 


antenna spacing. 
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Figure 28. Linear+circular array, SNR=15 dB, 7 elements, 2 
sources at (16°,15°), and (16°,20°) 
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Figure 29. Spiral array, SNR=15 dB, 7 elements, 2 sources 
at (16°,15°), and (16°,20°) 


oe) 


The circular array has a smoother spectrum and less 


spikes than linear+circular array, since the spacing of 


linear array part is larger than 0.54. From these 
Simulation results, we conclude that the larger antenna 
spacing will enhance angular resolution at the expense of 
more ripples on the spectrum. We can suppress the sidelobes 
and enhance resolution by increasing the number of elements 


in the x-y plane and increasing the radius of circle and 


then adding a linear part on z-axis with a spacing less 
than 0.52. 

Two elevation cases are considered. For ae large 
elevation angle (80 degrees), there are four. signals 
arriving at azimuth and elevation angles of (10°,80°), 
(Os cOs) eee oe 20 oO.) s with a  SNR=15 _ dB. The 


simulations are shown in Figure 34 for the circular array (7 


elements, =o and alvel Figure 35 iene the 


circular+linear array (3 elements, d=0.4A on z-axis). In 
both cases the sidelobes increase in the elevation 
direction, one Figure 35 shows lower sidelobes. From these 
results we conclude that the elevation angular resolution is 
worse when the elevation angle is large. This is a well 
known result that 1s associated with the increase in 
beamwidth that occurs with array scanning [Ref. 9]. 

For small elevation angles (5 degrees), there are four 
sources arriving at azimuth and elevation angles of 


(Gg = (00° (he ame 0,570, (270°, 5°) ead a, “SNR=15 sdb Ene 
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Simulations are shown in Figure 36 for the circular array (7 


elements, | ies es a and in Figure 27 for the 


circular+linear array (3 elements, d=0.44). Figure 37 
shows lower sidelobes in the plane of the sources. These 
sidelobes expand in the azimuth direction, implying that the 
azimuth angular resolution is worse when the elevation angle 
is small. This case is close to a broadside condition for 
the linear array and, conseguently, the resolution ability 
is essentially that of the linear array. Thus the resolving 
ability has been improved at low elevation by adding a 
linear array part (3 elements on the Z axis). 

Based on eRe results of the computer simulations and 
performance comparisons, we present some conclusions and 


recommendations in the next chapter. 
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Figure 30. Circular array, k,=1.52, SNR=15 dB, 7 elements, 2 
sources at (16°,15°), and (21°,15°) 
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Figure 31. lineart+circular array, d=1.54, SNR=15 dB, 7 elements, 
2 sources at (16°,15°), and (21°,15°) 
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Figure 32. Circular array, R,=1.5/4, SNR=15 dB, 7 elements, 2 
sources at (16°,15°), and (21°,15°), MUSIC method 
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Figure 33. Lineart+circular array, d=1.54, SNR=15 dB, 7 elements, 
2 sources at (16°,15°), and (21°,15°), MUSIC method 
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Figure 34. Circular array, SNR=15 dB, 7 elements, 4 sources at 
(10°,80°) , (90°,80°) , (180°,80°) , and (270°,80°) ,MUSIC method 





Figure 35. Linear+circular array, SNR=15 dB, 10 elements,4 sources 
at (10°,80°) , (90°,80°) , (180°,80°) , and (270°,80°) , MUSIC method 
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Figure 36. Circular array, SNR=15 dB, 7 elements, 4 sources at 
(10°,5°) , (90°, 5°) , (180°,5°), and (270°,5°), MUSIC method 
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Figure 37. Linear+circular array, SNR=15 dB, 10 elements, 4 
sources at(10°,5°) , (90°,5°) , (180°,5°) , (270°, 5°) , MUSIC method 
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VI. CONCLUSIONS 


Over the past few decades superresolution techniques 
have been developed to resolve multicomponent plane wave 
fields. The most noteworthy of these is the MUSIC 
algorithm. Others include PHD, Minimum norm, ESPRIT and 
conventional and Capon's beamforming. Each of the methods 
has its advantages and disadvantages with regard to hardware 
and computational aspects (array configurations and element 
spacing) and emitter parameters (number of emitters and 
their spatial distribution). 

This thesis has examined the performance of four of the 
estimation methods listed above for several linear, planar 
and volumetric arrays. Based on the formulations described 
in Chapter III, Matlab simulations of the methods were used 
to collect data for linear, circular, spiral and circular 
plus linear arrays. The emitter locations and spacings were 
varied, as was the signal-to-noise ratio. The data are 
presented as mesh surfaces, and the surface characteristics 
can be interpreted in terms of emitter resolution. 

For one-dimensional DF, the MRA provides the best 
dynamic range and resolution for a given number of elements 
and baseline. Since ESPRIT is generally applied to the 
spacing which is a known constant value, it is not directly 
applicable to MRAs. Using ESPRIT, the AOA can be obtained 


without a search technique, and in this regard is different 
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from the other four methods. Therefore, only the other four 
algorithms have been compared. 

The simulation resultS are in agreement with related 
two-dimensional DF performance simulations reported by 
Johnson and Miner [Ref. 10] for the MUSIC and Capon's MLM 
algorithms. The data in Chapter V leads to the following 
Goneclustonus- wiel) regard Go the circular “and weir,cular+tlinear 
a i Te 
(1) for high SNR levels, all four methods have about the 
same the performance, 

(2) elevation angular resolution is better than azimuth 
angular resolution for these three array configuration 
considered, 

(3) the azimuth angular resolution of the circular array is 
better than the spiral or circular+linear, 

(4) the elevation angular resolution of circular+linear 
array is better than spiral and circular+linear, 

(5) for multiple signal sources, the linear+circular array 
is better when antenna spacing is less than 0.5 wavelength, 
(6) larger element spacings result in better resolution, but 
there are more ripples and spikes; ambiguities will occur if 
the spacing becomes too large, 

(7) when the elevation angles of two signal sources are 
small, the azimuth angular resolution is worse than it is at 
large elevation angles, (because of slow falloff of the 


pattern in azimuth) 


SZ 


(8) when the elevation angles of two signal sources are 
large, the elevation angular resolution is worse than when 
elevation angles are small (because of slow falloff of the 
pattern in elevation). 

From the data it appears that the MUSIC algorithm is 
more stable, almost always correctly resolves individual 
emitters, and more importantly, has low sidelobes. However, 
it has been demonstrated that under some conditions MUSIC 
fails to provide the correct number of peaks [Ref. 11]. 

The circular array has best azimuth resolution but the 
results of linear+circular array are close to that of the 
circular array. Based on multiple signal sources’ the 


linear+cicular array is probably the better choice. When the 
antenna spacing is larger (say &k,=1.5 wavelength), the 


circular array is better. Having more elements on a larger 
diameter does provide some advantages. The distance between 
two adjacent antennas will be small to avoid high sidelobes 
and offers better azimuth resolution, especially for small 
elevation angle. 

Since we use subspace methods to figure out the AOA, 
the number of sources must be less than the number of array 
elements. Future work should consider how to deal with this 
limitation as well as improvements to the low signal-to- 
noise-ratio case. Any techniques that would allow a 
decrease in the search time for 2-D direction finding would 


also be valuable. 


63 


64 


APPENDIX A - THESIS MAIN PROGRAM 


clear all 

close 

Cv=version; 

it Strenumi cv 1a 5 


figure(l), 

bia ese jayne ao Ba 
(Suey le =. ex 4 
(Units 4s MOrmeligen Go, 


‘FPOsteron’, (0.0L -023 0.9 0.15 I, 

‘BackgroundColor 710.6 0.8 0.8], 

"HorizontalAlignment','left', 

'String', 'This program cannot be run on MATLAB 4.X,please use 
MATLAB 5.X'); 


pause(5), close,exit 
end 
load nsource 
msens=/; 
Snr=275- 


SnSOUrCe=2> 


ULCOU LT TOlL. Unies | points | : 
‘BaecrgrouncGolor ,;(O0xsw.6 (0.0), 
“VOMecol Zep L0y ? 
"Foreqgrounccolor’, (0 0 1); oo 
TPOSIELOn ,l2o.o 2oees 125. bios )|y 
‘Scring ' ; antenna number: ”, 

‘Steyle’ 7 Text’; ; 
‘Vag, octarltelext. 4 

Uicomecol | Unites’; pomts ;-%.s. 
UBackgrouncColor , [0.38 O28 U3], 
'POntoize , 10, 

PPOLecrOunccolLor 7 [140 07, ee 
MEGSLG TOM be 2. Oe oe to ee or ee 
YOLEING’ > Source Number: | 
SCYLe' 7* CORE", , 
"Tag, Stariclext |; 
UieonerOlLi Uribe 3) pOlme ~ro 
'BacCK@Groundgolo.*; |0<8 0-0 -0..6], 
VP OMpo bee’ 79) ce cc: 
“PeSit10On  ~ 122 208 (125-19 \, 
‘String’, signal to NemserRetictadsi:, 
ove a Lext : 
lag pp soetaclexe- a] 


msensl = wicentrol(*Units’;, points’, 


_ 


‘Backoroundcotor’ 0. L250, 

'‘Calipack* "ann: ; . 

"FOntSi ze’; 12, ; 
"POsteloe | lol gor JU. eons ols 
'Scering’ ,;ouUMZSeErlmsens) » 
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‘ScyLe’ - S@GisL | 
“"Eag' >; Beri ve es"); 
NScurcel = aecitrol( Unite . points’. 
“eckorouncColor’:, (0 1-0], 
‘Cal boaek 2. res... 
“EBOUNULSLZe pi es 
UVPeeiel On’; | l50.5.(2 ae 27 185i, 
"SEring'’,;numZstr(nseurce), 
Scyile"*, “eeue”,; 
NA Gioeeieme Pomc): jy 
Stirl = upcentrol( Units = peimts | 
‘BackgrounccColor (0° 1-0], 
‘Caliback’*;, arms; : 
YFOnt SL ze fia, 
“POSTE LOM » blo 0e: oo ne oll 
‘Suriung {~numZzstrtsnn), 
‘SOUVle | “EGrE") 23; 
leeone body Oiies ap oO rit a. ae 
'BacrqrouncCo lor {PU .o 0 .3"'0. 38 |, 
Perio ee eo 
LEOreqgroungcColorm oe like UeO), ; 
MEGSYVELON 3((40.9605 242209" 1002552 172-5798), 
‘SEPING = Sample poimnes: '; 
HOE bets Lex: je 
bnsamp=uicontrol('Units', ‘points’, 
WBaeckgrounceolor 70-1 O71. 
NG Sm aCie Gy seni fee eee 
LPOSiivonm. p~s(0502207 2422067926. 069 vl4 27759), 
PoEr LNG sg 4090. | 
roe y le eas, 
UIcOnETol( Units - poles, 
HEAiCKOrounaCe lor Lies 0 so Oe lip 
rEOmMpeolze alo, 
“hOre@g rounecoleor 4. (0) “Ori 2. 
"PESiiwon | [| bi eey oo 2672-379 54 16.756), 
‘SEring 7 Spacing. | 
Oey le’ > Lex” ),; 
bdspace= uicontrol('Units','points', 
(Beckoqrounecelor 7 |0 st s01, 
"Caollpack’, Inme. fsa 
'eOSr Lor, [252.150 268.071 26,009 1anZ ool, 
RSueimg 4 U4." y 
Poy be a Cale”); 
aw=clock; 


tname = ['Direction Finding Analysis---year:', num2str(aw(1)),’ 


MOonLh: 9) jj hUMZSEr(awl2 je, 7 —--cay'u »-MUMZSer (aw (ss) Ju 
set(gcf, 'name', tname) 


nsource=strZnum(get (nsourcel ,‘'string')); 
baspace!—“uleonerolt Units 72 poigts: 5 
rBackorouneecolor” 7 (Oulbsoi: 
Gal paek "Pama. "4. ccm: 
SPOSionom's([2e292 % 27 2], 
SET IG 4. O.4dy 
rouyle, "Sale ) 3 
bal= Uteonerol{ Unies’, pointes’; 
"Backgrounecolor’ 70 1:0]; 
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f 


Cent pack). inne 4s a 
"POSIULOn Gi gee ol eo2 12 
‘SerinG 5) Oy ; 
‘Style ) Sark )) 


MICONRLrOLT Unite | porte =, 
“BaGeagr7 oundcolor’;-(0ss 0.25020], 
"FOntS1 Ze 710, 
UPOoreqrounceolor | -( 0 Ori. 
"PONEWSI OnE ,. bola, 


PROSLE VOL pli lea ee eae le oes ei, 


rocring’, “l-D DIRECTION FINDPING =30 


PeeVle 4. Lext’, , 
‘lag’; Staticlexus jy 


MZCOntrolt Units | Doints: - 

‘Backgroundcolor’, [0.8 0.8 0,6, 

PFONLoVzZS 710, 

PEOUSOrDOUNICOLOT | ry. 
‘rOontWwelone bold’, 
PPOsteJol i 2Ov20-20 socal, poole, 
Poe rLne Gane les, 
~SOtyV le 3 Pest; , 
‘lag’, Statvetext4 *) > 


for n=1:8 

Dol=O.Conerol( nies @ polits |... 
"Backoroundtelor , (0.8 0238 0.3], 
"POS1t1ON (56.2571 64-56. 25). (n= 1) 
Pot etog 3. gone 
"Seyle’, text); 

end 


$clear bbl 
Sob1l=[] 
for n=l:nsource 
bbi (iy=uicontrol (* Units, pernts’; 
"BackeroundColor’, 10) 1 Oi, 
WCalibaeck” 7 inns” ; : 
MOStt ON filo 6.207 (ol oo. ol tne) 
we nai ele umme mal myc): om ad ears NAO as ag 
rOUyLe 2» COLe ') 3 
énd 
size(bbl); 
alee 
HiConoroi atts 2 polits 
'Baekoroundcoler', (0.6 028. 0.87, 
"FOntSize’ 7.9; 
"ForegroundColor','r', 
'FontWeight', ‘bold’, 


...90 deg', 


207 260.2 Jamo 


207 20.25 eee 


‘POStemen ~l2Z07o 1260.75 50,75 12.73), 


‘Sting’, Bo eangle:', 
rstyle, “Eext}; 


fOr 73-18 
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DebZin) = \UEcontrol( Units > “points 
"Baecxoroudgcolor’ , (02.6 0.6 0:3), : 
PGs eeiOn 71 o0.20F (90-62 2.25)" (ie-1)) 126275 26.25 Lae 
{Sereda ge pn tece 
‘Se yle- | 2ext oy 


ena 

clear bb2 

bp2z=[ J; 

for n=l:nsource 

Doz (in) = utcontrolL | *Unars; “pommes' 
"Backogrounacolor (0° io0), 
es er oye) cae el 9 Oe ho in a ae 
(PeOsitiron ; (ees2or (0-62.25) 9 (ial) 226 oe oo 
Scoring 7numz ser (10267 1), 
Seybe s COLe, ji: 

end 

ULCCNErO MN -WiatsS.~ > PoLlnTs. 


"BackgroundColor' -[0. 8 0. BS Ormisale 
UPOME SIZ 9 suey : 
BESO Ouro or, FOr Oe 
rPonmtweignt ' , “normal”, : 
Pe oS ei On.. alee One old 32. Zn ; 
String’ 7. 2—-)) DERECTION INDING: ELVangle.0 47250 mdeage “Ac angle: 
OIRnarane ts Cel = iain 
‘Style’ pile be 


Utcone sol Umita s. Pelmtse =. 
‘Bac grouncdeolor * hice 8 0. Sr Oe6i, 
'FontSize',9, 
“Pomewelon:', bola”, : 
‘Poa meron’, [2 #20G. 50.75 22.75), 
String» Bln angle.” 
fSty Le, “cexe); 

for n=1:8 

Dosim) = wWreenerol( Units. "polnes” 
‘Back@mounaeolor (0.8 OU. 85068 |, ; 
NEPOStsel On wu Sot (20, 75-05)" (i=l) LOSE 2or 2 oe Le 225), 
Sey ciauig (ey srs 
‘oEyl SC 5. Eexu 14 


end 

clear bb3 

bb3=[]; 

for n=l:nsource 

Coca -husconrcoln Unites points: 
So oe oe Pn 0 tag ee 8a Fe 
"Cal POaCkK sans se as 
"POSLELON j (ocean 7l. 75- 63) 24n=1) be) 26425 14.25) 
‘SEP LNG) pWUMZ seem 1)5 

‘Style *,. edit )3 
end 


bval= eo oe Bite S a DOLie Ss : 
‘BackgroundColor’ 10. 0 1), 
‘Colmigack 7. 1m"; 
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‘PONnGoLZe 4 LO, 

TFOREO rOuncdColor , (1-1 wots, 
PPOSLTELON 7g eo o2 eo) SOO 70d ore 4a) 
err Gate owe eo pe 2a ee ou ie 

'Style', 'popupmenu', 

'Tag', 'PopupMenul1'); 


clear bearingll 
clear bearingn 
clear ang 
msensll=strZnum(get( msensl, 'string')); 
snr=str2num(get(snrl ,'string')); 
dspace=str2num(get (bdspace, 'string')); 
nsamp=str2num(get (bnsamp, 'string')); 
thrv=get (bval, 'Value') 
if (thrv==1) 

thr=90; 
elseif (thrv==2) 

Chr=180; 
elseif (thrv==3) 

thr=270; 
elseif (thrv==4) 

thr=360; 
end 
for n=l:nsource 
bearingll(n)=str2num(get (bbl({(n),'string')) ; 


bearingn (n)=str2numiget (“bb2inj string: }); 

ang (m)=strzZnum(cee (bss (n) , String) )7 

end 

bearingll 

bearingn 

ang 

UPTcConerOol( UMmrEs 7 DOEMTS'. <x. 
"BackgroundColor, [0.757294 les02 7529410. 752040). 
"ButtonDownEcn', 'ctlpanel SelectMoveResize', 


"POstri0Om = 12zoo 151.5° 125-25. 64), 
‘olyile*, “Lrame-, 


‘lag’, Framel” ); 
MecOnt wel | Units "points, 
*"ButtonDownkcn', 'ctlpanel SelectMoveResize’, 


‘ot Paaeie., “Goan 5 
"FontSize',9, ae 
 POMeWwea rit ~ ~eaee "5 so: 
"POSPeron’ ; (Zoe ounce 24 2 oe loee Zoe 2045, 
"Strang ’ ~-<-Unitiorm and Non-unmilform Linear Array. 
“Lag PUStoureonw.: 17 
UL Contre. Units | worm: ss 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 
"Ca ltioaek , “Sprit”, 
'FOoRtoi1 ze 49, ha 
WE OMmEVictaite tO e, .-4% 
"POSUELVOM eee Deoe ele oe 
VString 4 “binear Array With BSOPRLE 
‘Vag so Pusnourtenm. ! ).; 
ULCOMEYOL( Units . pointes”, 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 
"Cal Vee". INvaEl” g 
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PEOuULo 2c 407, : 

'Foncwerghit', bold’, 

“E@elUVOn. ~icoe.0 166, 7/5711). 7 5255 i, 
"Sering ' , *MUSIC=resolution—= MOVIE 
Tag"; PUSADULEONS ©); 

PrcOnvnolL(- Units | points; 
'ButtonDownFcn', 'ctlpanel SelectMoveResize', 
"Callback', 'doag4', : 

'FoneAngle’,' italic’, 

EOE Nate... [IMEBS , 
ThOmMboUze 5 LU, 
‘Prone Wedoht 7 bold", ; 

PEOSt com si) te i oe) Oe oe ele 
‘String 4 2-D DIRECTION FINDING CIRCULAR 2RRAy 
hag go PUSshput eons 9 )> 

Liconvrolt “Units sy. ports 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 
"Cal Veaek’s) spira”™, : 

“POnecano le’ > rea lie. ; 
'PontName >." TIMES, 

MPOnt size 710; 
"Fonthwergnt’ , bola’, 

PEGS CLO mabe z <A eZ eo OZ le 
"String 5 27D DIRECTION FINDING SPIRAL ARRAY”, 
flag’, PusShbuttons )-; 

UPeone Lome Unrts | pores | 
‘ButtonDownFcn', 'ctlpanel SelectMoveResize', 

nGabipackit cil, 

Ponteangle, italic", 

"FontName*, ‘TIMES', 
SPOT Ote46 10. en 
TFOnEWeTOOL 4. DOL, 

TPOSTClOm) 71. b0Go 24a Zoe vO GZ Zod, 
'String','2-D DF LINEAR+CIRCULAR ARRAY', 
‘Tag’. Pusniburrons6. )% 

MEconerol| Units 7, pOines= ju. 
“Bae tor OundCOloGn +, (0.6 40 265 OU .oq-, 
PME@ne Size "4.9, 5% 
heme oroundcolor', (0 Oli, 
'POSiiEaen” (| Z2000759) 22.25.54 6.7586), 
"Swaing' , *LA=spacing: * , 
Psy ue *; text’); 

HLeEOn thom “Units | points | 
‘BackonroundColer', [Orcm0e cae. |e 

= “PontSi1ze "9, : 
VPOreoroundcolor 3,[0 707 I), ; 
VOCS ir Om 740 Ow Oo Oo acon 4 a Ow OC, 
Strang bA-ank iNo. 
Style >. texte 4); 


ULeonenol( Unies |. pointes =| 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 
'Gallback"; close”, 
Pontsaze y iZ, 
Ome Vic oie 7015 Od Cc a: 
MBGie tne dG eal ee ee eee ees | 
“SECU 7 (Co LOos | 
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‘Tag’, Pushbutton? “1; 


Miceontrol( Units - polis | 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 
Vea ge lock » 


(POMetolze 7412, 
"FontWeight', 'bold', 
PROsieLon 7, (29. Oe Zoot oe Oe toe ole 
Pee tic ae COC. 
Tae pe UshioUecon.) 


MiSOntrOl(s Vitis . DOLGhS jr. 
'BackorounoColor’., (0.6 0.3 0.3). 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 

'FontName', 'times', 


TORE Ze pa, 
"rormeuveighit’;, “bold”, 
"POreorounacolor , | 1-007, ’ 
“POSPELOR Yilze0 o2o 27s 5 leo oe 
'String', ‘Antenna Geometry', 
'Style'’, 'text', : 
‘Tag 4 Steet 2cCl1esto ); 
MEGOntrol( = Units; BoIneEs 7. 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 
"Callback','sp2', 
"FOntsi2ze '. 9, 
"FontWeight', 'bold', 
“POS UUUON 2100) 24 oole ene ee oS 


"String, Nen-Uunitorm linear array; 
“tag 4 PusSnouErone, 1 
Uvegnorolt "Unites, pOines:. | 


"ButtonDownFcn', 'ctlpanel SelectMoveResize', 
‘Callback’, sel’; 

‘Fonts. Ze >. 9; 
'FontWeight', 'bold', 

PPGSit hor , 1290.25 221.75 1229 25e22. 5, 


‘String’, "Uniform Linear Array', 
Teg. Pushbuctons <j. 
URCOncrol( y Unites. jy Colmes. 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 


MC cinLOgek 2 Oe optus 
*FontSize’, 10, 
‘FontcWeregnt.pemte 7 . 4. 
POSTEO pile foe ol ee es oe ey 
‘String » 2-D DF Brray ; 
“Tag. Pushbuttons "); 
Plecterol ( UaTes S oolnes 
"ButtonDownFcn', 'ctlpanel SelectMoveResize', 
‘Callback ) Nelpwitm((sce, cel) | 
"FOR GoOLZe yo, 
‘roncwelgne , bOla”. 45% 
"POSitacm', (346.5) 2504/5 247.29 4712 
MOLE Ing 7 BELP; 
‘Taq 7 -UsShbuelons-. || 


ULGOnE rol ¢ UmiLs ; “poincs,; 
TBackoreuna olor’, (0.8. 0.8 0. 2], 
'ButtonDownFcn', 'ctlpanel SelectMoveResize', 
'FontName', 'times', 


yg 


'Foncervzel7 13, 
'Fontweight’,' bold’, 
‘PeoregreuncColer:, [1 0 0], ave 
‘Posie Lon; |209.5 Ale 2a 2 or ese oe 
"String , Antenna Pattern ; 
'Sayle’; Rext is 
MEControl Ones; Dottie s 4 
'ButtonDownFcn', 'ctlpanel SelectMoveResize', 
"Cal loaeck rs, “celap’; 
“WOmeS ize >. 5 ae 
"hongwWweigne . oceld"' ,; 
PEPOed Gl On pike 9 we hee oe oy 
YSEring » Laneartcircilar array). 
ULEONEr OL Unavs. , points | 
'ButtonDownFecn', 'ctlpanel SelectMoveResize', 
‘Callback ~ Claw; 
'SOntolze toy 7 
bPontwWwetohe "74: bold>, 
UP Os Me lOlls dy wy co 2 oie et ae oy Geto 


USE Lingo eireularvariray &)> 
str = {' This program was created by LIN, KU-TING, TAIWAN' 
' Date:9/98 


Thesis eopie 

' Comparsion of superresolution algorithms with different array 
geometry’ 

Ww fOr redio. Gireceion Ldn | 


Thesis adviser:Prof. David C.Jenn' 
" second reader:Prof. Pace Phillip’ 


i ' 


We evaluate the effectiveness of various array designs and 
processing' 
'" method for two-dimensional direction finding and resolving 
MmulLcreple* 
" emitters. Using five different spatial spectrum methods to 
predict ehe’ 
" direction of arrival (DOA) ' 
Discussion Of Topics: * 
‘ 1. spatial spectrum estimation methods' 
) MUSIC(multiple signal classification' 
J eESeRET.: 
) PHD(Pisarenko harmonic decomposition) ' 
) minmum-norm' 
) Capon beamformer ' 
rray Geometry' 
) Unweeonm Linear array (i—v).” 
) non-uniform linear array(1-D) ' 
) “Circular osarray (2-0)? 
) spiral array(2-D)' 
Noliicamterrcular array (z2-D)" 
3. Antenna Pattern' 
(a) circular array’ 
(5). limneartcireular array’ 
4. Angular Resolution' 


a2 


3 


‘HOW TO USE." 


1. Input parameter' 
(a)antenna number:less than 10 for non-uniform linear 


array' 

(b) source number:less than antenna number' 

: (Ce) SNR(GB):signal to noise ratio,” 

(d)sample points: the more points,the more accurate 
results! 

(e)spacing:1-D DF antenna spacing(normalized to wave- 
length) ' 


2-D DF the distance from center point to 
element (X-Y plane) ' 

(f{) LA spacing:antenna spacing for linear array part (Z- 
axis)of linear+circular array' 

: (g) LA-ant No:antenna number for linear array part(Z- 
axis)of linear+circular array' 

(h)1-D direction finding angle range:0~180 deg' 

(1)2-D DF EL angle:0~90 degree' 

; (j)2-D DF AZ angle:0~360 degree' 

note:choose the AZ searching range --> push popup 


menu! ' 

: only four choices:90,180,270,360 degree. 

; 2. When you change the "source number",please wait one 
moment ! ' 

: 3. Check antenna array geometry and antenna pattern' 

: 4. Choose 2-D or 2-D DF --> Push button,sceeiwiar you ger 

SS ee een 

1. The 2-D DFs take some time to search the AOAs,' 

please be patient to wait,the word "HOLD ON!!" will show 
lena 

; 2. ESPRIT method just for uniform linear array:20 runs' 

: 3. MUSIC-resolution-movie DEMO --> four tagerts fixed, one 
moves. ' 

sample points=200' 

: 4. This program CANNOT be run on Matlab 4.X, at least 
Matlab 5.X ' 

: 5. The default antenna number of lineart+circular array is 
EO (=743)i 

; 6. The larger radius of circular array,the better 
resolution’ 

7. When the elevation angle of AOA is small(5 degree), the 
sidelobe' 

/ expand on the azimuth direction.’ 

8. When the elevation abgle of AOA is large(80 degree),the 
sidelobe’ 


: expand on the elevation direction.' 

: 9. When the radius of circular array is larger(1.5 
wavelength),there' 

are more spikes, supressing the spikes by adding linear 
array on the Z-axis' 

“ote =" DOA’; 
enelpwin(str, ttl); 
teset(O265, 1.a0Z2fsetstr(108),  fontname’, * symbol", “fomusize 7.20) 
axis (Ole }) 
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APPENDIX B - SPECTRUM CALCULATION FOR CIRCULAR ARRAY 


format compact 
bearing=bearingn; 
msens=msens1ll1; 
msource=nsource; 


def msens = 8; 
dere msouUrce = 2, 
smat = zeros(nsamp,msource); for -n = 1:msource 
Uo= rend (hsamp, 1) —C.5; 
SMac ls; 1) = —siagn(u) =~ oq l=2-aosi(u))= 
end 
%  ¢€+------------ > 
bearing = bearing * pi / 180; S degrees to rad 


phi=2*pi/msens; 

qi=[0:msens-1] *phi; 

for nnn=1: msource 
Ciratnin,. —qi; 

end 

for nb=l:msens 
cb(:,nb)=bearing'; 


end 
amat2=((2*pi*dspace*cos(cira-cb)))'; 
for ns=l:msource 
a3(:,nsS)=amat2(:,ns) *Sin(ang(ns) *pi/180) ; 
end 


amat 3=exp(j%*a3) ; 


te) 


% the noisefree sensor signal 
smat2=smat*amat3.'; 


[nsamp, msens] = size(smat2) ; 
nsource=msource; 
$snr=input ('snr='); 


namp=1/ (sqrt (2)*10*(snr/20) ); 
for n=l:msens . 


nmat(:,n) =randn(nsamp,1)+j*randn(nsamp, 1); 

end 

S$nmat = nmat - ones(nsamp,1) * mean(nmat); 
smat2=smat2+namp*nmat; 

ymat=smat2; 

=== estimate the spatial cross-cumulant matrix ---- 
cmat = zeros(msens,msens); 

ymat = ymat - ones(nsamp,1) * mean(ymat); 6 remove mean 
cmat = conj{ymat' *~ ymat) / nsamp; 6 correlation matrix 


a Getermine Number cl “SCurces ==——s52=—-------.--— 


fumat, smat, vmat] = svd(cmat); 
svec = diag(smat); 


Formac. COmMpact 


° 


= estimate Dearing Speci ra. SS sas ea y eS 9 == SS 56s eo eee a ee 


prmat = zeros(msens,msens) ; 
for i=l:msens 
Lf (svyecii) > 0) 


prmat = prmat + vmat(:,i) * vmat(:,i)' / svec(i; 
end 
end 
prvec = prmat(:,1); 
dtheta=1; 


theta = [O:dtheta:thr]' * (pi/180); 
ve=(021:901]~ (617180) ; 
pp=(length(theta)); 
for nnl=l1:pp 
Ciralinni;. }—aqi; 
end 


for nbl=l:msens 
ebl-(:, nb1l)=theta: 
end 


omega = 2*pi*dspace*cos(ciral-cbl); ¢<-=-<---=--------= 
BRAAARR HHH HHARRRAAAA 


mth = lengenteneta); 

mic = zeros(mth,1); 

qvec —.umac(l, be asource)-..5 

gmat = umat(2:msens, 1:nsource) ; 

arvecc= ()lo5° = -cenj tomate) 7 qvec yo (1 = gvec’*gvec) j7 

tmat = vmat(:,nsourcet+l:msens) ; 

$wvmat = tmat * diag(ones(msens-nsource,1) ./ svec(nsource+l:msens)) * 
tmat' ; 

Ete: = ehaca “ema. ? 


nv=length(va); 
if thr==360 
for bb=l:nv 
for J} =l men 


steer = expt] * omegqaik,:) ~sanm(vai(bb)) 1: 
Smusl(k, bb) = lafwmstcer Mie senate Estee: = 
end 
end 
else 


for bb=l:nv 
£O2 kK —lemeh 


steer = 23D" omega tk, =) ssi (va (bb) ) 7; 
beaml(k,bb) = 1./abs(steer' *inv( cmat)* steer); 
sminl(k,bb) = 1./ abs(arvec.' * steer)%*2; 
smusl(k,bb) = 1./ (steer' * tmat * steer); 
spisl(k,bb) = 1./ abs(steer' * vmat(:,msens) })%2; 
end 
end 


end 
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Val=|=90212 90] “i417 130)7 
mtl=length(val); 


eae display the angular spectra ----- 


theta = theta*180/pi; 
Va=Va~ 13/1; 
val=val’*180/pi; 


top=0.95; spacing=0.02; 
left=0.80; 
btnWid=0.1; 
btnHt=0.04; 
bottom=0.05; 


figure (1) 
DLeoneroll Units 4 polnt Ss. | 
"BackgroundColor, (0.3 023-028); 
ROnCoLZe yaks, . 
TPOrecTrOoungce lor 40.6. Veo (ea |, 
YPOSit1On 4270 200.100 1914 
Potrane 4. HOlD. CN 
TS ye go tem 
‘Tag’, 'StatieTextz*)]; 
fa=ger(0, 'curmentiireure’); 
figure(fat+l), 


as=10*loglO(abs( (smusl+eps)')); 
asl=10*loglO(abs( (smuslt+eps)')); 
smesh( theta',va,10*loglOQ(abs( (smusl+eps)'))); 


bb=10*1log10(abs((smusl+eps)')); 
bias=min (min(bb) ); 
as=max (max (bb) ); 


if thr==360 

mesh( theta',va,bias+10*logl0Q(abs( (smusl+eps)'))); 
Stitie({ MUSIC, Cireular erray,;SNR=",nuMmzZser(snr,3), «Yedceus-¢— 
Spacing.” ,numzser (dspace, 3) |)> 

Zlabel('power(dB)'); 

xlabel('AZ(deg)') 

ylabel('EL(deg) ') 


axist((0 2O"thry 0 90-20 asi); 
$xlabel({'ANTENNA NUMBER=',int2str(msens)]); 
Sylabel(['source number=',int2str(nsource)]); 
zlabel ('power') 


break 

else 

as2=10*logl0(abs( (spislt+eps)')); 
as3=10~logl0O(abs( (sminl+eps)')); 


as4=10*logl0(abs( (beaml+eps)')); 
subplot (221),mesh( theta',va,bias+10*loglO(abs( (smuslt+eps)'))); 
zlabel ('power(dB)'); 
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tetio( MUSIC*)>xlabel "Az(deq) *) 
ylabel('EL(deg) ') 


axis([0O 90*thrv 0 90 -20 as]); 
oxlabel({'ANTENNA NUMBER=', int2str(msens)]); 
soylabel(['source number=',intZ2str(nsource)]); 
szlabel ('power (dB) '') 


bb=10*logl0(abs((Sspislt+eps)')); 
bias=min(min(bb)); 
as=max(max(bb));> 


Subploti(Z222),mesh © theta’, va; biastlO*logl0(abs (spisl+eps) 72))4-- 
behets PDs. 

zlabel ('power(dB)'); 

xlabel ('AZ (deg) ') 

ylabel ('EL(deg) ') 


axis([0 90*thrv 0 90 -20 as}); 
$xlabel ('Az-angle'"); 


bb=10*logl0(abs((sminlteps)')); 
bias=min(min(bb) ); 
as=max(max(bb)); 


subplot (223),mesh({ theta'’,va,bias+10*loglO(abs{ (smini+eps)'))); 
Eeelet Mint =norm a 

Zlabel ('power(dB)'); 

xlabel ('AZ(deg) ') 

ylabel ('EL(deg) ') 


axis([0 90*thrv 0 90 -20 as]); 
Sylabel ('El-angle'); 
bb=10*logl0 (abs ((beaml+eps)')); 
bias=min(min(bb)); 
S$as=max(max(bb)); 


subplot (224),mesh( theta',va,abs (bias)+10*logl0(abs((beamlt+eps)'))); 
title('Capon MLM’) 
Zlabel ('power (dB) ' 
xlabel ('AZ(deg) ' 
ylabel ('EL(deg) ') 
end 


Vs 


AZ=bearing*180/pi; 
EL=ang; 


Stitle(['ANTENNA NUMBER=',int2str(msens),' SOURCES 

NUMBER=', int2str (msource) ] ) 

$xlabel(['AZ=',int2str(AZ(li:msource)),'' degree']}) 

svilabel{ |" Bh="Ganezstrien(l:msource)),- degree']) 

tname = ['CIRCULAR ANTENNA ARRAY,SNR=',int2str(snr)]; 
set (gcf, 'name', tname) 
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APPENDIX C - SPECTRUM CALCULATION FOR SPIRAL ARRAY 


format compact 

bearing=bearingn; 

def msens = 8; 

Geil msource = 2; 

msource=nsource; 

msens=msensll; 

smat = zeros(nsamp,msource); 
for n = l:msource 
o> rand (usamp, 1) -0. 5; 
smat(?,;n)= “siqniu) =* Lo¢ (i=2ap61u)); end 
6 the propagation matrix 
bearing = bearing * pi / 180; 


phi=2*pi/msens; 
ql=(O:msens-1] *phi; 
dsp=[l1:msens] *dspace* (1/msens) ; 
Sesp—(1 2 5 11°13 16°) ie 

for nnn=1l: msource 
Clr (inn, :)=qi? 

alta, +) =asp; 

end 

fer nb=lsmsens 
cb(:,nb)=bearing'; 


end 
amat2e=((2"*pi"*Ci. cos (Cira—-ch))) 
for ns=l:msource 
a3(i,ns)=amatZ(:,ns)*sin(ang(ns) *pi7180) ; 
end 


amat 3=exp(j*a3); 
smat2=smat*amat3.'; 
[nsamp, msens] = size(smat2); nsource=msource; 


namp=1/ (sqrt (2) "10° (snr/20)); 
for n=l:msens 
nmmat(:,n) =randn(nsamp,1)+j*randn(nsamp,1); 


end 


smat2=smat2+namp*nmat; 


ymat=smat2;%----- estimate the spatial matrix ---- 
cmat = zeros(msens,msens) ; 
ymat = ymat - ones(nsamp,1l) * mean(ymat); 3% remove mean 
cmat = conj(ymat' * ymat) / nsamp; 
fumat, smat, vmat] = svd(cmat); 
svec = diag(smat); 
format compact 
+ estimate bearing spectra -------- 
prmat = zeros(msens,msens) ; 


for i=l:msens 
Lt Svec{(i) > 0) 
prmat = prmat + vmat(:,1) * vmat(:,1)' / svec(i); 
end 
end 
prvec = prmeat(:,1); 
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dtheta=1; 
theta = [O:dtheta:thr]' * (pi/180); 
=U 21290} ea LSU); 
pp=(length(theta)); 
for nnl=l:pp 
Giral(nnd, = os; 
d(nnl,:)=dsp; 
end 
for nbl=l:msens 
cb1(:,nbl)=theta; 


end 
Omega -="2 "pia. * Cosi(el ral-ebl); 
mth = length(theta); 
mic = zeros(mth,1); 
gvec = umat(l,l:nsource).'; 
gmat = umat(2:msens,1:nsource) ; 
arvec = [1 ; - conj(gmat) * gvec / (1 ~- gvec'*gvec) 
tmat = vmat(:,nsourcetl:msens); 
wvmat = tmat * diag(ones(msens-nsource,1l) ./ 
svec(nsourcet+l:msens)) * tmat'; 
tmat = tmat * tmat'; 
nv=length(va); 
ae Seoecancieanig Le .———————— S 
bOr. bb=leny 
One he 
steer = @€xD (7) 7 “Omegatk, =) °"*sinivaloo) ) }- 
beamit(k, bo)-=14/--abs (steer* *inv( cmat).~* steer) ; 
smirk, bbje= 1.7 vabstarvec.  ““steer) 72. 
smusl(k, bo) = 127 (steer * tmatz — steer); 
Sei sik, 55) 1./ abs(steer' * vmat(:,msens) )%*2; 
end 
end 


welts POO) toi 1 B07 

mtl=length(val); 

jae Orsplaystne angular spectra. ----— 
theta = theta*180/pi; 

va=va*180/pi; 

val=val*180/pi; 


BOp=0 957; spacing—U .02; 
left=0.80; 
btnwWid=0.1; 
btnHt=0.04; 
bottom=0.05; 
figure(l) 
ULGOReroli mars. peincs”;, . « 
Tuacmenoundcolon regio ge! 0. 8 0. Cale, 
[RONG orce. wo, «42 
Tice Celie LOE OP ce Oe oa ae 
PPOSLE1ON pal Ly orezode VOOe Irs |, 
VS Peer 7 HO WDe eM hs 
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ie 


OL Le a tex ey . 
"Tag, statliclexiZ 


fa=get (0, 'currentfigure') 
figure(fat+l), 
tname = ['SPIRAL CIRCULAR ANTENNA ARRAY,SNR=',int2str(snr)!; 
set (gcf, 'name', tname) 
asl=10*loglO(abs( (smuslteps)')) 
as2=10*loglO(abs( (sSpislt+eps)')); 
as3=10*loglO(abs( (sminlt+eps)')); 
as4=10*loglO(abs( (beaml+eps) ')) 


f 
1 


, 


bb=10*log10(abs((smusl+eps) ')); 
bias=min(min(bb)); 
as=max (max (bb) ); 


if thr==360 


mesh( theta',va,bias+l0O*loglO(abs( (smuslt+eps)'))); 
sti tle( (MUSIC, Spiral circular arrayr~jSNR= ,nuUmMZStr (snr, 3) ,* ,max- 
Spacing: ",numZzstr (dspace; 3) |); 


axis([0O 90*thrv 0 90 -20 as]j); 
Zlabel ('power(dB)'); 

xlabel ('AZ(deg) ') 

ylabel ('EL(deg) ') 


else 

subplot (221),mesh( theta',va,bias+10*logl0(abs( (smusl+eps)'))); 
$title('*music'); 

ax2zst[O SO*thry 0 90) =20 asl); 


bb=10*1log10(abs((spislt+eps) ')); 

bias=min(min(bb)); 

as=max (max (bb) ); 

zlabel ('power(dB)'); 

ticle MUertC’ )fxlabelt Azideg) *) 

ylabel ('EL(deg)‘) 

subplot (222),mesh( theta',va,bias+10*loglO(abs( (spislt+eps)'))); 
Lace (PHD. 

zlabel ('power(dB)');xlabel('AZ(deg)');ylabel('EL(deg)') 

axis (tO SO enry-& 90 =2U as); 


bb=10*logl0(abs((sminlt+eps)')); 

bias=min (min (bb) ); 

as=max(max(bb)); 

S$xlabel ('Az-angle'); 

subplot (223),mesh( theta',va,bias+10*loglO(abs( (sminlteps)'))); 
title('Mini-norm') ; 

zlabel ('power(dB)'); 

title('MUSIC');xlabel ('AZ(deg)') 

ylabel ('EL(deg) ') 


axis([O 90*thrv 0 90 -20 as]); 

$xlabel(['uniform circular array,SNR=',intZstr(snr) )); 
sylabel('’El=angle’); 

bb=10*1logl0(abs((beamlt+eps)')); 

bias=min(min(bb)); 
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subplot (z2.4),mesh{ theta',va,abs (bias) +10*1l0g10 (abs ((beamlt+eps!')}); 
title('Capon MLM'),; 

zlabel('power(dB)'); 

xlabel('AZ(deg) ') 

ylabel ('EL(deg) ") 

end 

AZ=bearing*180/pi; 

EL=ang; 


mtitle(['SPIRAL ANTENNA NUMBER=',int2str(msens),' SOURCES 
NUMBER=', int2str(msource) ] ) 
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APPENDIX D - SPECTRUM CALCULATION FOR CIRCULAR+LINEAR ARRAY 


dspacel=str2num(get (bdspacel, 'string')) 
im=str2num(get (bdl, 'string')) 
bearing=bearingn; 
msource=2; 
if nsource>6 
msource=6; 


else 

msource=nsource; 
end 

smat = zeros(nsamp,msource) ; 

for n = l:msource 

u = rand(nsamp,1) -0.5; 
Ssmat\:,n)— —Sign(u) =a* Log(1l—=-2*abc{u) )}> 

end 


° 


S$ the propagation matrix 
bearing = bearing * pi / 180; % degrees to radian 


msensm=msensll; 
phi=2*pi/msensm; 


ql=[O:msensm-1]*phi; 
Cira=ones (msource,1)*ql; cb=bearing'*ones(1,msensm) ; 


amat2=((2*pi*dspace*cos(cira-cb)))'; 
amacz2=({(2*o1* (02 lm-1) **dseacel)) 7 


fOr *ns—) 2msource 


a3(:,ns)=amat2(:,ns)*sin(ang(ns)*pi/180) ; 
a32(:,ns)=amat22*cos (ang(ns)*pi/180); end 
a3=[a3;a32]; 


amat3=exp(j*a3); 


oo 


the noisefree sensor signal 
smat2=smat*amat3.'; 


oe 


--=a00 white Nolse=———. 
[nsamp, msens] = size(smat2); 
nsource=msource; 


namp=1/ (sqrt (2)*10*(snr/20) ); 
for n=l:msens 


nmat(:,n) =randn(nsamp,1)+j)*radn(nsamp,1); 
end 
nmat = nmat - ones(nsamp,1) * mean(nmat); 


smat2=smat2+namp*nmat; 
ymat=smat2; 


2) 


Noein estimate the spatial cross-cumulant matrix ---- 
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cmat = zeros (msens,msens); 


ymat = ymat - ones(nsamp,1: * mean(ymat); _ remove mean 


cmat = conj(ymat' * ymat) / nsamp; 2 correlation matrix 


a. aes determine MAimber Ol Sources \=--===-——-5--_. | 


[umat, smat, vmat] = svd(cmat); 
svec = diag(smat); 

s=[1 2 3]; 

dtheta=l1; 


format compact 


'S- = =—=-- I ee ee ee ee 


prmate= Zeros (msens,msens); 
for i=l:msens 
if (svec(i) > 0) 


Prva = prmat + VMer(s 2) *-vmeet 11 2 sovecia): 
end 
end 
wEevec = pumati: 7); 
theta = m0 sdrneta rin 4" (pr/ls0, 


Va=[(0: 1790] * (61/180); 
pp=(length (theta) ); 
eGireali=ones (ppl) al, 
cbl=theta*ones (1,msensm) ; 


omega = 2*pi*dspace*cos(ciral-cbl); %<---@@@e@e~** 

mth = length(theta); 

mic = zeros(mth,1); 

Gvec = tae () .linscurce).~ ; 

gmat = umat(2:msens,1:nsource) ; 

arvec = [1 °fe= con7z (omat) * qvec / “(1 = qvec'*gvec) |]; 

tmat = vmat(:,nsource+l:msens); 

wvmat = tmat * diag(ones(msens-nsource,1) ./ svec(nsource+1:msens) ) 
tmat  ; 

tmat = tmat * tmat'; 


nv=length(va); 
if thr==360 
for bb=l:nv 
ror kK =lemth 


steerl = exp(j * omega(k,:)'*sin(va(bb))); 

Steen7— Gx *“Ccosiva(b>)) ~2*e1*dspacei~ (O:im-1)"); 
steer=[steerl;steer2]; smusl(k,bb) = 1./ (steer' * tmat 
end 

end 
else 


for bb=l nv 
for k =l:mth 


steerl = exp(j * omega(k,:)'*sin(va(bb))); 
Steen7— exp) CoOsiva (b> )y) "Zi “dspace!” (O:tm-1)]*) 7 
steer=[steerl;steer2]; 
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* steer); 


* 


beaml(k,bb! = 1./abs(steer' *inv( cmat)* steer); 
Sminl(k,bb) = 1.7/7 abs(arvec.' * steer)*2; 
smusl(k,bb) = 1./ (steer' * tmat ~* steer); 
spisl(k, bb) 1./ abs(steer' * vmat(:,msens) )%2; 
end 

end 

end 


Val=(=903 1" 50) * (oy leo) 


mtl=length(val); 


= display the angular spectra ----- 
theta = theta*180/pi; 


va=va*180/pi; 
val=val*180/pi; 


top=0.95; spacing=0.02; 
left=0.80; 
btnWid=0.1; 
btnHt=0.04; 
bottom=0.05; 
figure(l) 
Umeonreol(” Units | bormts; 
reackerouncdcColor |)( 0.8.0 0c]. 
(Ponte o Ize 1) o- : 
'PorneqroundColor', (0.3. 0.8_0.3], 
‘Position’, (178.264 100 19), 
rotring’ » HOLD ONT", 
'Style', 'text', 
“Tag; oSbatilclext2 7; 


fa=get (0, 'currentfigure'); 
ProuLre (tearl); 


asl=10*loglO(abs( (smuslteps)')); 
bb=10*1logl0(abs((smusl+eps)')); 
bias=min(min(bb) ); 
as=max (max (bb)); 


if thr==360 

mesh( theta',va,bias+10*logl0(abs( (smusl+eps)'))); 

statle( (“MUSIC, Linéearrcitewiar array,oNR=—",oUNZSEr isnt, 3) 7 .e- 
Spacing: ',numZstr (dspace, 3),',l-space>”* ,numZstr (dspacel, ss), 7 L-eant 
No: . fumZstrisms) js 

axis([0O0 90*thrv 0 90 -20 as]); 

$xlabel(['ANTENNA NUMBER=',int2str(msens) ]); 

Sylabel(['source number=',int2str(nsource) }); 
S$zlabel('power(dB)'); 

zlabel('power(dB)'); 

xlabel ('AZ (deg) ') 


SS 


ylabel/'EL(deg)') 


else 

as2=10*loglO(abs( (spisl+eps)');); 

as3=10*logl0O(abs( (sminit+eps)')); 

as4=10*logl0(abs( (beaml+eps)')); 

subplot (221),mesh( theta',va,bias+10*loglO(abs# ‘smuslteps) "I§}; 


axis(({0O 90*thrv 0 90 -20 as]); 
zlabel('power(dB)'); 

title ('MUSIC');xlabel ('AZ(deg)') 
ylabel ('EL(deg) ') 


bb= Pe eee cerns reek) 2 
bias=min (min (bb) ); 
as=max (max (bb) ); 


subplot (222),mesh( theta',va,bias+10*loglO(abs( (spislteps)'))); 
Zlabel ('power(dB)'); 

ertie( PHD’) 7;xlabel( Agideg)) 

ylabel ('EL(deg)') 


axis(({0 90*thrv 0 90 -20 as]); 
bb=107 loolO(abs (t(sminiteps) “)); 
bias=min(min(bb)); 
as=max (max (bb) ); 


Subp louetz2s), mesh tnete 7 va, orast1 0" logl0 (abst (sminitens) 9?) )e 
title, *Mani= ee oe 
zlabel ('power(dB)'); 


ee ties ) 
ylabel ('EL(deg) ') 


avs ( [0 S0*thry 0 90 =20>4s5])- 
tname=([{'linear+circular,SNR=',int2str(snr)]); 
Sylabel ('El (degree) '); 


bb=10*logl10 (abs ( (beaml+eps)')); 
bias=min(min(bb)); 
$as=max(max(bb)); 


subplot (224),mesh( theta',va,abs(bias)+10*logl10(abs((beamlt+eps)'))); 
tvele<" Capon MLM); 
Zlabel ('power (dB) ' 
xlabel ('AZ(deg) ') 
ylabel('EL(deg)') 
end 
AZ=bearing*180/pi; 
EL=ang; 
tname=([{'linear+circular, SNR=',int2str(snr)]); 
set(gcf, 'name',tname) 


); 
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APPENDIX E - SPECTRUM CALCULATION FOR LINEAR ARRAY 


clear nmat 
[nsamp, msens}] = size(smatl); 
nsource=msource; 


namp=1/ (sqrt (2)*10*(snr/20)); 
for n=l:msens 
nmat(:,n) = randn(nsamp,1)+j*randn(nsamp, 1); 


end 
smati=smatl+namp*nmat,; 


ymat=smatl; 


s-——-_ estimate the spatial cross-cumulant matrix ---- 


cmat = zeros(msens,msens) ; 
ymat = ymat - ones(nsamp,1l) * mean(ymat); 6 remove mean 
le determine mumbper Ol SOUurces: —=——==—===—==—---==— 
cmat = conj(ymat' * ymat) / nsamp; 6 correlation matrix 
foo > CeLermitic MUMDer Of ‘SCUYCES Soha assse>—4-5-se— 
[umatc, Smear, vnet | = Ssyvyd(cmeat i: 
svec = diag(smat); 
ae i estimate: Dearing SPeeia, 2A ee ee 
prmat = zeros(msens,msens); 


for i=l:msens 
at wiSswec (1) > 70) : 
prmat = prmat + vmat(:,1) * vmat(:,i)" / Svec(i)-; 
end 
end 
prvec = prmat(:, 1) - dtheta—O0.5; 


theta = [-90:dtheta:90]' * (pi/180); 

omega = 2*pi*dspace*sin (theta) ; 

mth = length(theta); 

mic = zeros(mth,1); smus = mlc; spis = mlc; seig = mlc; par = 
nic; 

beam = mle, Smin = mic; 


Qvec = UMae(!l, banseource) .”; 


gmat = umat(2:msens, 1:nsource) ; 
arvec = [1 ; - conj(gmat) * gvec / (1 - gvec'*gvec) j; 
tmat = vmat(:,nsourcet+tl:msens) ; 


oy 


wyvmat = tmat * diagiones(msens-nsource,1) ./ svec(nsource+l:msens) ) 


tmat'; 
tmat = tmat * tmat'’; 


for k =l:mth 


steer = exp(j * omega(k) * [O:msens-1]'); 
obeam(k) = abs(steer' * cmat * steer); 
beam(k) = 1./abs(steer' * inv(cmat) * steer); 
smin(k) = 1./ abs(arvec.' * steer)%2; 
smus(k) = 1./ (steer' * tmat * steer); 
spis(k) = 1./ abs(steer' * vmat(:,msens) )%2; 
end 
ee ee display the angular spectra ----- 


theta = theta’180/pi; 


spec = abs( [seig, smus, spis, mlc, par, smin, beam] ); 
6-==c >>> normalize to abs-max of unity for display only 
spmax = max(spec); 

spmax = ones(l1,7) ./ spmax; 

spec = spec * diag(spmax); 

figure (2) 


S$subplot (421), stem(svec), title(tname) 

SObplOE (271) pOplorteneta, SO" logLOtspec( sy 3) ) 7s) 
title, Unitorm Jinear array) 

ylabel ('power (dB) ') 

Geext (1 PHD | 

subpler(423), plot leheta, 0 logit (spect 2) jr) 
ylabel ('power(dB)'), gtext('MUSIC') 

SUbDLeGEd42>)) plete neua,, 10 leg) Utspec(=. 6) 47m) 
ylabel ('power(dB)'),gtext('Mini-norm') 

Sucve or (427) 7 plou(tneca,. 0" loglotspeci(?, 7.5, =), nold 


ylabel ('power(dB)'),gtext('Capon MLM'),xlabel ('AOA(degree) ') 


tname = ['UNIFORM ANTENNA and NON-UNIFORM ANTENNA LINEAR ARRAY' J 


set(gcf, ‘name’, tname) 
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APPENDIX F - SPECTRUM CALCULATION FOR NON-UNIFORM LINEAR 


ARRAY 
clear nmat 
[nsamp, msens] = size(smat2); 
nsource=msource; 
$snr=input('snr='); namp=1/(sqrt(2)*10*(snr/20)); for n=l:msens 
nmat(:,n) =randn(nsamp,1)+7*randn(nsamp,1); 
end 
smat2=smatZ2+namp*nmat; 
ymat=smat2; 
=== estimate the spatial cross-cumulant matrix ---- 
cmat = zeros(msens,msens); 
ymat = ymat - ones(nsamp,1l) * mean(ymat); 6 remove mean 


cmat = conj(ymat' * ymat) / nsamp; 


correlation Metrix 


o> = GOeternine Mimier Ol SOULESS Soe — a= =— 25 -- oe 5 


fumat, smat, vmat] = svd(cmat); 
svec = diag(smat); 
NolG Ore 
$-—------- estimate: bearing. specera ==—-—-——-=—-—=-=---—---____ -_-_ + 
prmat = zeros(msens,msens) ; 
for i=l:msens 


aE (Svecti1)>> 0) 


prmat = prmat + -Vmeat yt) = wmac(s,4) 7 / “svectiis 
end 
end 
prvec = prmat(:,1); 
dtheta=0.5; 
tmeta = (|[=SGedtheta:90)' *~ (pa71s0); 
omega = 2*pi*dspace”~ sin (theta); 
mth = length(theta); 
mc = zeros(mth,1); 
smus = mlc; spis = mlc;seig = mlc; par = mic; 
beam = mlc; smin = mic; 
gvec = UMmat (1, Limsource). 77 
gmat = umat (2:msens,1:nsource) ; 
arvec = [1 ; - conj(gmat) * gvec / (1 - gvec'*gvec) ]; 
tmat = vmat(:,nsourcet+l:msens) ; 
wvmat = tmat * diag(ones(msens-nsource,1) ./ svec(nsource+tl:msens)) * 
tmat'; 
tmat = tmat * tmat'; 


for. =i mck 
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steer = exp(sqrt(-1) * omega(k) * asp'); 


wpeam(k}) = abs(steer' * cmat * steer); 
beam(k) =l./ abs(steer' *inv( cmat) * steer); 
smin(k) @ 1./ abs(arvec.' * steer)%*2; 

smus(k) = 1./ (steer' * tmat <* steer); 
spis({(k) = 1./ abs(steer' * vmat(:,msens) )%2; 


a eas one display the angular spectra ----- 


theta = theta*180/pi; 
spec = abs( [smus, spis, smin, beam] ); 


i —--_— normalize to abs-max of unity for display only ----- 
spmax = max(spec); 


spmax = ones(1,4) ./ spmax; 
spec = spec * diag(spmax); 
EECQuULed 2) 


Subp lOu(4 22) 7 tolor(theta, VOrlogl lO spects4 2) 05 
SSeb GCay DOSEt LOM al C400? 02/5. 0.23905'0 217 |e 
title('Non~-uniform linear array'),ylabel('power(dB)'),gtext('PHD'); 


subplot (424), 

ploutthetas 10" hegll(spec(:7 1) )) pylabeh{ power (as) orexti MUSIC"). 
SUbppLOE (426) 7 plot (ineta,. l0*logl0<seec( 4,33 ))> 

ylabel ('power (dB) '),gtext('Mini-norm' ) 

SSCUCCCae position ~ (0206 "WU s25 0e35.07 17 ia 

S$title('Minimum norm method') 


SUDDPOR(428) ,7p loc theta, -W0* locl0ispects,4))); 
ylabel ('power(dB)'),xlabel('AOA(degree)'),gtext('Capon MLM’) 
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